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1  EXECUTIVE  SUMMARY 

Our  work  during  this  contract  period  began  with  emphasis  on  electromagnetically  induced 
transparency,  slow  light,  and  its  application  to  making  ultra- long  biphotons.  As  the  program 
evolved,  we  increased  emphasis  on  novel  techniques  to  produce,  measure,  and  utilize  both 
long  and  short  biphotons. 

We  begin  by  summarizing  the  key  property  of  time-energy  entangled  biphotons.  This  is: 
if  an  observer  at  point  A  chooses  to  measure  the  frequency  of  an  arriving  photon  he  will  then 
know  to  high  accuracy  the  frequency  of  the  photon  which  will  be  measured  by  an  observer 
at  point  B.  But  instead,  if  the  observer  at  point  A  chooses  to  measure  the  time  of  arrival  of 
a  photon  at  his  location,  he  will  then  know,  again  to  high  accuracy,  the  time  of  arrival  of 
the  photon  at  point  B.  The  accuracy  of  these  measurements  is  not  limited  by  the  Heisenberg 
uncertainty  principle. 

We  turn  next  to  what  is  meant  by  long  and  short.  Typical  biphotons  as  generated  by 
nonlinear  optical  crystals  have  temporal  lengths  in  the  range  of  between  0.1  ps  and  10  ps. 
Such  photons  are  not  resolvable  by  presently  existing  photo  detectors.  These  detectors, 
measure  whether  or  not  a  photon  is  present,  but  may  not  be  used  to  examine  the  functional 
form  of  the  photonic  wave  packet. 

By  using  the  techniques  of  slow  light  and  working  under  AFOSR-ARO  and  DARPA 
support,  in  2005  our  group  demonstrated  the  first  method  of  generating  temporally  long 
biphotons.  The  length  of  these  photons  is  controlled  by  the  group  delay  in  the  nonlinear 
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media  and  in  early  experiments  resulted  in  photons  with  a  length  of  about  40  ns,  A  key 
accomplishment  of  the  present  program  was  the  extension  of  this  length  to  photons  whose 
length  could  be  continuously  varied  from  50  to  900  ns.  Of  importance,  the  line  width  of 
these  photons  was  less  than  the  natural  linewidth  of  the  rubidium  vapor  that  was  used 
to  produce  them.  This  is  important  because  optical  nonlinearities  when  produced  using 
electromagnetically  induced  transparency  continue  to  increase  in  the  subnatural  linewidth 
regime  and  in  the  absence  of  dephasing  may  be  made  arbitrarily  large.  This  should  soon 
allow  the  demonstration  of  nonlinear  optics  with  single  photons. 

We  mention  a  surprise  that  occurred  during  this  work.  In  the  course  of  observing  long 
biphotons  with  a  length  determined  by  the  slow  optical  group  velocities,  we  found  that 
the  photonic  wave  packets  had  a  sharp  leading  edge  spike  on  their  front  edge.  Following 
a  suggestion  by  Dan  Gauthier  we  recognized  that  this  front  edge  spike  is  a  Sommerfeld- 
Brillouin  precursor.  This  observation  is  important  because  it  clarifies,  for  both  slow  light 
and  fast  light,  that  information  will  always  be  transmitted  at  the  speed  of  light  in  vacuum. 

In  January  2008  we  recognized  that  we  had  the  capability  to  modulate  single  photons  for 
the  first  time.  To  do  this  we  used  the  Stokes  photon  of  a  biphoton  pair  to  set  the  time  origin 
for  electro-optic  modulation  of  the  wave  function  of  the  anti-stokes  photon.  With  the  time 
origin  determined,  the  modulator  could  arbitrarily  modulate  either  the  amplitude  or  phase 
of  the  anti-stokes  photon.  The  technique  therefore  provides  the  technology  for  studying  the 
response  of  atoms  to  shaped  single-photon  waveforms  on  a  time  scale  comparable  to  the 
natural  linewidth  of  target  atoms. 

The  next  step  in  our  work  in  modulating  biphotons  was  the  development  of  a  method  for 
measuring  their  length  using  slow  detectors.  The  essential  idea  is  that  modulation  in  the  time 
domain  followed  by  slow  integration  constitutes  a  Fourier  transformation.  The  experimental 
technique  is  to  measure  the  coincidence  count  rate  between  single  photon  counting  modules  as 
a  function  of  an  applied  sinusoidal  frequency.  The  inverse  Fourier  transform  of  the  data  then 
yields  the  biphoton  waveform.  Though  this  experiment  was  a  proof  of  principle  experiment, 
ultimately  it  could  be  used  to  measure  wavepacket  profile  of  biphotons  when  sufficiently  fast 
photo  detectors  are  not  available. 
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We  are  now  working  on  an  important  extension  of  our  work  on  modulating  single  pho¬ 
tons.  This  is  the  application  of  spread  spectrum  techniques  at  the  single  photon  level.  Spread 
spectrum  is  well  known  in  the  communications  industry  as  a  technique  for  avoiding  inter¬ 
ference  and  jamming,  and  at  times  increasing  information  capacity.  Our  work  is  the  first 
demonstration  of  this  technique  to  single  photons. 

We  turn  next  to  our  work  on  short  biphotons.  Our  first  contribution  was  a  PRL  entitled 
Chirp  and  Compress  :  Toward  Single-Cycle  Biphotons.  In  this  work  we  described  a  method 
for  generating  time-energy  entangled  photons  with  a  spectral  width  that  exceeds  an  octave, 
and  of  compressing  their  spectrum  to  produce  biphotons  whose  temporal  length  is  a  single 
optical  cycle.  The  elements  of  the  technique  are  the  suggestion  for  using  parametric  down 
conversion  in  a  periodically  poled  material  to  spontaneously  generate  pairs  of  entangled 
photons  whose  instantaneous  frequencies  are  chirped  in  opposite  directions  and  the  use  of 
the  non-local  nature  of  entangled  photons  to  allow  the  dispersion,  as  experienced  by  one 
photon,  to  cancel  out  the  dispersion  of  the  second  photon  and  to  compress  the  biphoton 
wave  packet. 

Our  first  experimental  work  in  this  area  was  the  development  of  a  novel  resonant  sum 
frequency  generation  technique  (PRL  February  2009).  This  work  demonstrates  a  rather 
amazing  effect  where  we  take  two  single  photons  that  each  have  a  broad  spectral  linewidth 
and  sum  them.  We  find  that  though  each  photon  has  a  broad  linewidth,  the  linewidth  of 
the  sum  frequency  photon  is  as  narrow  as  that  of  the  pumping  laser.  Also,  the  output  (sum) 
power  is  linear  rather  than  quadratic  in  input  power.  These  effects  both  derive  from  the 
quantum  behavior  of  single  photons. 

In  the  course  of  studying  nonlocal  dispersion  compensation  as  described  above,  we  recog¬ 
nized  that  there  should  be  a  new  quantum  effect  that  we  have  termed  as  nonlocal  modulation. 
Assume  that  single  and  idler  photons  pass  through  sinusoidal  phase  modulators  located  at 
different  locations.  These  modulators  are  driven  at  the  same  modulation  frequency  and  are 
connected  by  cable  such  that  their  relative  phase  may  be  varied.  After  passing  through 
the  modulators  the  single  and  idler  photons  are  dispersed,  for  example  by  a  prism,  and 
the  relative  positions  of  the  single  and  idler  photons  are  correlated.  We  find  :  When  the 
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modulators  are  run  with  the  same  phase  the  modulation  depths  add;  when  they  are  run  in 
phase  opposition  the  modulation  depths  subtract.  Two  distant  modulators  with  the  same 
modulation  depth  and  opposite  phase  therefore  have  the  same  frequency  correlation  as  when 
both  modulators  are  absent.  This  effect  is  entirely  quantum  mechanical.  Mathematically 
it  results  because,  quantum  mechanically,  one  adds  probability  amplitudes  before  squaring, 
while  classically  one  squares  before  adding. 

As  this  contract  ends,  we  are  beginning  work  in  the  area  on  nonlinear  and  quantum  optics 
at  x-ray  wavelengths.  To  a  great  extent  we  are  motivated  by  the  new  1.5  Angstrom  free 
electron  laser  that  is  now  operating  at  Stanford.  We  anticipate  experiments  to  demonstrate 
frequency  doubling  of  this  laser  with  the  further  objective  of  using  SHG  as  a  correlator  and 
diagnostic  for  the  laser  itself. 


2  Subnatural  Linewidth  Biphoton  Generation  with  2-D 
MOT 

The  most  significant  experimental  advance  in  cold  atom  EIT  research  project  during  perfor¬ 
mance  period  is  that  we  observed  subnatural  linewidth  biphotons  with  our  MOT  generation 
II,  which  is  our  new  platform  for  biphoton  generation  in  a  series  of  modulation  studies  at 
single  photon  level.  Figure  1  shows  both  experimental  configuration  (geometry)  and  mech¬ 
anism  of  parametric  paired-photon  generation.  In  contrast  to  previous  MOT  constructed  in 
early  grant  period  (W911NF-04-1-0105),  the  new  2-D  ^^Rb  MOT  has  a  cigar  shaped  atom 
cloud(~1.7  cm  long  and  an  aspect  ratio  of  25)  and  consequently  a  large  optical  depth  in  the 
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(b) 


|4>5P3^(F=3) 


Figure  1:  Upper:  Biphoton  generation  in  a  double-A  system,  (a)  Experimental  configura¬ 
tion.  Fi  and  F2  are  narrow-band  optical  frequency  filters,  (b)  ®^Rb  energy  level  diagram. 
In  the  presence  of  counter-propagating  pump  {cUp)  and  coupling  (lUc)  beams,  Stokes  {tOg  )and 
anti-Stokes,  (tUas)  photons  are  generated  into  opposing  single-mode  fibers.  Lower:  2-D  MOT 
apparatus.  The  vacuum  cell  is  6  cm  size  ceramic  structured  octagon.  The  cell  is  located  in 
the  middle  of  water  cooled  trapping  coil(racetrack-shaped  cage). 

longitudinal  direction;  moreover,  its  cylindrical  quadrupole  trapping  field  results  in  minimal 
longitudinal  magnetic  field  gradient  and  hence  greatly  reduces  the  inhomogeneous  Zeeman 
broadening  of  the  m-states  of  the  55i/2  level.  The  experimental  cycle  comprises  4.5  ms  of 
trapping  time  and  0.5  ms  paired  photon  generation  window.  At  the  end  of  the  trapping 
cycle,  the  rubidium  cloud  is  prepared  in  55i/2  level  by  turning  off  the  repumping  laser  0.3 
ms  before  turning  off  the  trapping  laser;  counter-propagating,  circularly  polarized,  cw  pump 
{(Vp)  and  coupling  (uJc)  lasers  are  subsequently  turned  on  and  phase-matched,  paired  Stokes 
(cVs)  and  anti-Stokes  (was)  photons  are  spontaneously  generated  and  propagate  in  opposite 
directions  as  shown  in  the  figure. 
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Figure  2:  Biphoton  wave  packet  data  for  three  slow  group  delay  cases.  The  upper  row  plots 
are  anti-Stokes  EIT  scan  data(o)  and  EIT  fit(blue  curve),  the  lower  row  plots  are  paired- 
photon  coincidence  count  data(-l-)  and  predicted  wave  packet  shape  with  time  bin  width  of  1 
ns  for  800  s(lower  row).  Propagation  delay  of  anti-Stokes  pulse(red  traces  in  inserts  of  lower 
row  plots)  are  also  presented  for  three  slow  group  delay((Tg)  cases.  Experimental  parameters 
are:  Left:  (r^  ~  50ns)  OD=7,  flc  =  4.2O713,  flp  =  I.I6713,  and  Ap  =  48. 67713.  Middle: 
{Tg  ~  320ns)  OD=53,  tic  =  4.2O713,  f2p  =  I.I6713,  and  Ap  =  48.67713.  Right:  (r^  ~  900ns) 
OD=53,  a  =  2.35713,  ftp  =  I.I6713,  and  Ap  =  48.677i3. 


The  optical  depth  of  the  2-D  ^^Rb  MOT  can  be  varied  up  to  62,  which  gives  us  enough 
parameter  space  to  verify  the  relation  between  the  optical  group  delay  and  the  length  of 
the  biphoton  waveform.  Figure  2  shows  sets  of  anti-Stokes  EIT  scan  and  paired-photon 
coincidence  counts  for  three  anti-Stokes  EIT  group  delay  cases,  which  is  controlled  by  varying 
the  optical  depth  and  the  coupling  laser  power(rp  ~  (27i3/|r^cp)A^(TL,  where  NaL  is  the 
optical  depth  and  713  is  the  dephasing  rate  of  level  |3)).  The  trapping  laser  used  for  these 
experimental  runs  has  a  power  of  160  mW,  a  beam  diameter  of  2  cm,  and  is  red  detuned  by 
20  MHz  from  the  |55i/2,F  =  3)  ^  |5P3/2,F  =  4)  transition.  A  repumping  laser  is  locked 
to  the  |55i/2,F’  =  2)  15^3/2,7^  =  2)  transition,  has  a  power  of  80  mW,  and  overlaps  one 

of  six  trapping  beams.  The  pump  laser  is  circularly  polarized  (c"),  has  a  1/e^  diameter  of 
1.46  mm,  and  is  blue  detuned  from  the  |1)  — >  |4)  transition  by  146  MHz,  i.e.  Ap  =  48. 67713. 
The  coupling  laser  is  circularly  polarized  (cr+),  has  a  1/e^  beam  diameter  of  1.63  mm  and 
is  on  resonance  with  the  |2)  — >  |3)  transition.  The  counter-propagating  pump  and  coupling 


7 


beams  are  collinear  and  set  at  a  2  degree  angle  from  the  longitudinal  axis  of  the  MOT.  The 
Stokes  (a”)  and  anti-Stokes  ((7“^)  photons  are  coupled  into  opposing  single  mode  fibers  after 
passage  through  A/4  wave  plates  and  polarization  beam  splitters  (PBS).  The  Stokes  and 
anti-Stokes  fiber  coupling  efficiency  is  70%  and  the  1/e^  waist  diameter  of  their  foci  is  220 
/xm. 


Group  delay  (ns)  Optical  depth  (OD) 

Figure  3:  (a)  Measured  correlation  time  vs  measured  anti-Stokes  group  delay.  The  solid  line 
is  a  linear  least  squares  fit.  (b)  Paired  counts  in  a  1  ns  bin  in  800  s  as  a  function  of  the 
optical  depth. 

As  shown  in  Fig.  2,  the  temporal  length  of  biphoton  wavepacket  generated  in  high  OD  case 
are  much  longer(by  more  than  one  order)  than  previously  reported  in  our  phase  I  research. 
It  also  means  that  the  generated  biphoton  has  much  narrowed  bandwidth.  The  predicted 
biphoton  packet  waveforms  in  lower  row  plots  of  Fig.  2  are  computed  with  all  parameters 
obtained  from  the  EIT  measurements  and  vertically  scaled  to  fit  the  experimental  data.  The 
calculated  biphoton  linewidths  are  9.66,  2.36  and  0.75  MHz  respectively.  These  linewidths 
are  comparable  to  the  measured  EIT  bandwidths  and  in  the  latter  two  cases  are  less  than 
the  6  MHz  natural  linewidth  of  Rb  D  line.  Thanks  to  the  2-D  MOT,  which  is  designed  to 
have  increase  in  OD  and  decreases  in  dephasing,  the  biphoton  generation  now  is  in  the  linear 
group  delay  regime  where  Tg  >  Tr  and  the  correlation  width  directly  follows  group  delay  time 
as  shown  in  data  plot  of  Fig.  3(a). 

Having  taken  into  account  the  filter  and  etalon  transmissions,  the  fiber  to  fiber  coupling 
efficiency,  the  detector  quantum  efficiencies  and  the  duty  cycle,  for  the  conditions  of  Fig.  2, 
we  observe  a  total  of  3213,  31674,  and  22000  paired  counts  in  800  seconds,  which  correspond 
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to  generation  rates  of  1275,  12569,  and  8730  pair/s,  respectively.  Higher  generation  rates 
can  be  achieved  by  increasing  the  pump  laser  power.  With  Qp  =  6.88713,  the  paired  photon 
generation  rates  are  4.0  x  10^  and  9.0  x  10^  pair/s  at  OD=17  and  30,  respectively.  Figure  3(b) 
shows  that  the  number  of  paired  counts  varies  linearly  with  the  optical  depth.  Though 
the  generation  rate  per  spectral  bandwidth  varies  as  the  square  of  the  optical  depth,  the 
bandwidth  reduces  linearly  with  this  depth  leading  to  the  linear  dependence.  Experimentally, 
under  optimum  conditions,  we  observe  74%  of  the  Stokes  photons  to  be  paired.  We  also 
observe  that  all  of  the  correlation  data  violate  the  Cauchy-Schwarz  inequality  by  as  much  as 
a  factor  of  11600. 

It  is  noticeable  that  there  is  a  sharp  peak  at  the  leading  edge  of  the  correlation  data 
generated  for  high  optical  depth  case  as  shown  in  lower  plots  in  middle  and  right  column  of 
Fig.  2.  This  feature  is  Sommerfeld-Brillouin  precursor  type  as  in  the  case  for  propagating 
classical  wave  packets.  The  physical  picture  is  the  following:  the  detectors  register  biphoton 
coincidence  counts  versus  the  time  t  =  tas  —  tg]  the  earliest  portion  of  the  biphoton  wave 
packet  comes  from  the  high  frequency  portion  of  the  spectrum,  which  is  not  in  the  range  of 
large  group  delay,  and  is  thought  of  as  the  Sommerfeld  precursor.  At  slightly  later  times  the 
low  frequency  Brillouin  components  arrive  at  the  detector  and  beat  with  the  simultaneously 
arriving  high  frequency  components.  Though  precursors  are  now  understood  in  the  optical 
region  and  have  even  been  observed  long  ago  with  correlated  gamma-ray  photons,  our  work 
reports  the  first  observation  of  precursors  as  measured  by  single  photon  correlation. 

The  2D  MOT  we  have  developed  is  likely  to  allow  immediate  improvements  in  single 
photon  read-write  techniques,  in  EIT-based  quantum  memory,  and  in  nonlinear  optics  with 
cold  atoms.  Each  of  these  areas  requires  the  same  ingredients  as  demonstrated  here;  i.e.,  low 
dephasing  rate,  subnatural  linewidth,  and  high  optical  depth.  For  example,  the  efficiency  of  a 
nonlinear  process  continues  to  improve  as  the  photon  linewidth  is  reduced  below  the  natural 
linewidth,  and  is  ultimately  limited  by  the  dephasing  rate  of  the  non-allowed  transition. 

For  more  details  please  see:  1)  Shengwang  Du,  Pavel  Kolchin,  Chinmay  Belthangady, 
G.  Y.  Yin,  and  S.  E.  Harris,  ’’Subnatural  Linewidth  Biphotons  with  Controllable  Temporal 
Length,”  Phys.  Rev.  Lett.  100,  183603  (May,  2008);  2)  Shengwang  Du,  Pavel  Kolchin, 
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Chinmay  Belthangady,  G.  Y.  Yin,  and  S.  E.  Harris,  ’’Observation  of  Optical  Precursors  at 
the  Biphoton  Level,”  Optics  Letters  33,  2149  (September,  2008). 


3  Theory  of  EIT-  based  Paired  Photon  Generation 

Our  theoretical  model  is  highlighted  on  paired  photon  generation  in  the  double-lambda 
atomic  system.  We  first  reexamine  the  conditions  required  for  the  system  to  operate  in 
the  group  delay  regime.  We  find  that  the  optical  depth  of  the  atomic  system  has  to  be 
sufficiently  high  in  order  to  avoid  the  filtering  of  the  generation  bandwidth  of  paired  photons 
by  the  transparency  window.  Second,  we  have  extended  the  theoretical  treatment  beyond 
the  ground  state  approximation.  This  allows  us  to  properly  include  and  analyze  the  effect  of 
Langevin  noise  fluctuation  on  the  atomic  system  and  solve  the  problem  of  its  return  to  the 
ground  state  after  the  emission  of  Stokes  and  anti-Stokes  photon  pairs.  We  also  addressed 
the  important  questions:  1)  '‘What  are  the  effects  of  the  optical  thickness  of  the  atomic 
sample  on  paired  and  single  photon  generation?”  2)  “Does  every  Stokes  photon  have  its 
paired  anti-Stokes  photon?” 

With  low  parametric  gain  and  high  optical  depth  we  show  that  the  system  can  produce 
highly  correlated  photon  pairs.  The  shape  of  the  intensity  correlation  function  and  the 
emission  bandwidth  depend  on  the  coupling  laser  Rabi  frequency  and  the  optical  depth  of  the 
atomic  sample.  Compared  to  SPDC,  paired  photon  generation  in  the  double-lambda  atomic 
system  is  affected  by  Raman  gain  in  the  Stokes  channel  and  EIT  in  the  anti-Stokes  channel. 
EIT,  through  the  absorption  at  the  poles,  cuts  the  emission  bandwidth.  In  order  to  enter 
a  regime  where  the  EIT  window  is  sufficiently  large  and  therefore  the  emission  bandwidth 
is  controlled  to  a  large  extent  by  the  phase-matching  process  in  the  presence  of  large  group 
delay,  the  optical  depth  of  the  atomic  sample  has  to  be  larger  than  10.  High  optical  depth 
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substantially  reduces  the  influence  of  Langevin  noise  fluctuations  and  Raman  scattering  on 
paired  photon  generation  so  that  the  Stokes  and  anti-Stokes  photons  are  generated  mostly  in 
pairs.  We  therefore  suggest  the  use  of  a  cigar  shaped  atomic  cloud  with  high  optical  depth 
in  the  longitudinal  direction. 

For  more  details  please  see:  P.  Kolchin,“Electromagnetically-lnduced- Transparency-Based 
Paired  Photon  Generation,”  Phys.  Rev.  A  75,  033814  (March  2007).  Kolchin’s  dissertation 
chapter  4  “Theory  of  EIT  based  Paired  Photon  Generation”  is  included  in  this  report  as 
Appendix  C. 


4  Electro-Optic  Modulation  of  Single  Photons 

We  have  demonstrated  conditional  shaping  of  single  photon  waveforms  using  electro-optic 
modulators.  We  use  the  Stokes  photon  of  a  biphoton  pair  to  set  the  time  origin  for  electro¬ 
optic  modulation  of  the  wave  function  of  the  anti-Stokes  photon.  This  technique  allows 
arbitrary  control  of  both  phase  and  amplitude  of  single  photons.  A  key  requirement  for  this 
modulation  scheme  is  that  the  temporal  length  of  biphoton  wavepacket,  which  is  represented 
by  correlation  time  in  experiment,  is  much  longer  than  modulation  response  time  designed 
for  experiment.  The  subnatural  line  width  biphoton  source  described  in  the  previous  section 
producing  biphotons  with  correlation  times  adjustable  in  the  50-900  ns  range  is  ideally  suited 
for  this  application.  Our  single  photon  counting  modules(about  350  ps  of  timing  resolution) 
and  data  system  has  fast  enough  temporal  resolution  to  observe  AM  modulation  in  this 
work. 

Figure  4  shows  the  schematic  of  the  experiment.  We  use  counterpropagating  cw  pump 
and  coupling  lasers  to  generate  time-energy  entangled  pairs  of  Stokes  and  anti-Stokes  photons 
which  propagate  in  opposite  directions  and  are  collected  into  single  mode  fibers  as  described 
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Figure  4:  Schematic  of  paired  photon  generation  and  conditional  modulation.  A  Stokes 
photon  detected  by  an  SPCM  sets  the  time  origin  for  shaping  the  anti-Stokes  photon  with 
an  electro-optic  modulator.  To  within  the  accuracy  of  the  SPCM,  this  allows  shaping  of 
both  the  amplitude  and  phase  of  the  anti-Stokes  photon. 

in  section  A.  The  detection  of  the  Stokes  photon  at  D1  sets  the  time  origin  for  firing  the 
function  generator  that  drives  the  electro-optic  modulator  which  in  turn  modulates  the  anti- 
Stokes  photon.  Verification  of  the  single  photon  nature  of  the  modulated  anti-Stokes  photon 

/ty\ 

is  done  using  a  50-50  beam  splitter  and  detectors  D2  and  D3.  Their  coincidence  data  G23TT) 
is  expected  to  show  a  dip  at  r  =  0  proving  nonclassical  nature  of  modulated  anti-Stokes 
signal. 

The  electro-optic  amplitude  modulator  consists  of  phase  modulators  in  both  arms  of  a 
Mach-Zehnder  (MZ)  interferometer.  The  degree  of  phase  control  in  both  arms  depends  on 
the  type  of  the  modulator.  We  use  a  z-cut  modulator  that  requires  =  1-75  volts  to 
cause  the  tt  phase  shift  required  to  go  from  minimum  to  maximum  transmission  and  can 
be  operated  at  a  maximum  frequency  of  10  GHz.  One  port  of  the  output  beam  splitter 
of  the  MZ  interferometer  is  terminated  so  that  the  portion  of  the  photon  wave  function 
that  is  not  transmitted  is  lost.  In  general,  if  a  Stokes  photon  is  detected  at  time  ti,  and 
the  modulator  is  activated  conditioned  on  this  detection  then,  in  the  Heisenberg  picture, 
the  anti-Stokes  operator  at  the  output  of  the  modulator  is  related  to  the  input  operator 
by  &out{h)  =  f  9it2,t2)ain{t2)dt2.  If  there  are  no  dispersive  elements,  then  to  within  an 
unimportant  phase  factor  we  may  write  df^atih)  =  '^{^)dinit2)  •  The  correlation  function  in 
the  presence  of  the  modulator  is  related  to  that  in  the  absence  of  the  modulator  by 

Gm\^)  =  (4.1) 
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With  the  biphoton  wavefunction  given  by  +  r)  the  modulated  (conditional)  single 

photon  wavefunction  is  +  r).  We  adjust  the  bias  voltage  at  the  input  of  mod¬ 

ulator  so  that  the  output  of  the  modulator  m(r)  is  related  to  the  input  voltage  V{t)  by 
m(r)  =  sin  [(/>(r)]  exp  [za(/>(r)],  where  (/)(r)  =  7rV(r)/(214),  and  a  is  a  phase  modulation 
parameter.  For  a  z-cut  amplitude  modulator  as  used  here  a  =  0.75,  but  may  be  eliminated 
by  using  an  x-cut  modulator. 
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Figure  5:  D1-D2  coincidence  counts  in  a  1  ns  bin  as  a  function  of  the  delay  between  Stokes  and 
anti-Stokes  photons,  (a)  Modulated  (A)  and  unmodulated  (□)  waveforms,  (b)  Waveforms 
with  Gaussian  (□)  and  rising  exponential  (a)  shapes.  The  experimental  data  (□,  A)  were 
collected  over  2000  s.  The  solid  curves  for  cases  (a)  and  (b)  are  plotted  from  theory.  The 
inset  in  part  (a)  is  the  scope  trace  of  the  output  voltage  of  the  function  generator. 


The  principal  experimental  results  of  this  work  are  shown  in  Figure  5.  In  part  (a),  the 
modulation  signal,  shown  as  an  inset,  is  a  set  of  two  square  pulses.  Of  importance,  there 
is  no  vertical  scaling  between  the  modulated  and  non- modulated  waveforms.  In  Fig.  5(b), 
we  show  photons  modulated  with  two  different  waveforms.  In  the  first  case  the  modulator 
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is  driven  with  a  gaussian  pulse.  In  the  second  case  we  design  the  function  generator  output 
so  as  to  compensate  for  the  nonlinear  distortion  of  the  modulator  in  such  a  way  that  the 
output  of  the  modulator  is  an  exact  rising  exponential. 

We  define  the  retrie\^l  efficiency,  Sr  of  a  paired  photon  source  as  the  probability  to  gener¬ 
ate  a  single  anti-Stokes  photon  on  th('  condition  that  its  paired  Stokes  photon  is  detected.  For 
the  non-modulated  photon  we  measure  Sr  =  3.5%.  When  losses  at  the  beamsplitter,  modu¬ 
lator,  filters,  fiber  to  fiber  coupling  and  detector  efficiency  are  backed  out,  this  corresponds  to 
a  retrieval  efficiency  of  55%.  For  the  modulated  photons  the  measured  retrieval  efficiencies  of 
the  two  square  pulses,  the  rising  exponential  and  gaussian  waveforms  are  Sr  =  1.3%,  0.61% 
and  0.9%  respectively.  With  losses  backed  out  these  efficiencies  are  21%,  9.4%  and  11.2% 
respectively. 


Figure  6:  (color  online)  Conditional  three-fold  correlation  function  gcond(^)  ^  ^  function 
of  the  Stokes  rate  for  unmodulated  (□)  and  modulated  (A)  single-photon  generation.  The 
dashed  curve  shows  the  theoretical  limit  for  in  the  absence  of  excess  light  scattering 

(see  text). 

Since  single  photons  incident  on  a  beamsplitter  must  go  into  one  output  port  or  the  other, 
in  the  ideal  case  where  there  are  no  two-photon  events  and  there  is  no  excess  scattered  light, 
we  would  expect  no  three-fold  coincidences  at  the  detectors.  A  measure  of  the  quality  of 
heralded  single  photons  that  quantifies  suppression  of  two  photon  events  is  given  by  the 
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conditional  correlation  function: 


(2)  ^r\\  "^123 -^1  /A 

^cond(O)  -  JVi2A^i3  •  ^  ^ 

Here  A^i  is  the  number  of  the  Stokes  counts  at  Di,  N12,  and  N13  are  the  number  of  two-fold 
coincidence  counts  within  a  time  window  Tc  at  detectors  Di,  D2  and  Di,  D^;  and  N123  is 
the  number  of  three-fold  coincidence  counts  within  this  same  time  window. 

In  Figure  6  triangles  and  squares  show  measured  versus  Stokas  rate  with  and 

without  modulation.  The  modulation  is  done  with  the  same  signal  as  in  Fig.  5(a).  We  set 
Tc  equal  to  the  nominal  length  of  the  unmodulated  biphoton  (285  ns).  At  a  Stokes  rate  of 
2.2  X  10^  sec“^  which  corresponds  to  Qp  =  0.26  Fa,  we  obtain  5cmd(0)  ~  0-2  i  0.04  and 
ffcondW  ~  0.21  ±0.07  for  the  unmodulated  and  modulated  waveforms  respectively.  The  fact 
that  the  measured  is  less  than  0.5,  (the  limiting  value  for  a  two  photon  Fock  state), 

is  indicative  of  the  near-single  photon  character  of  the  light  source. 

Because  there  is  a  small  probability  for  the  parametric  down  conversion  process  to  generate 
multiple  pairs  of  biphotons,  even  in  the  absence  of  spurious  light  scattering,  the  conditional 
correlation  function  is  not  zero.  The  dashed  curve  In  Fig.  6  shows  the  theoretical  prediction 
for  the  conditional  correlation  function  that  results  from  such  multiple  scattering  events. 
Because  of  light  scattering  from  both  the  pump  and  coupling  lasers,  the  experimental  curves 
lie  above  this  limiting  value. 

We  perform  two  control  experiments:  In  the  first  we  remove  the  30  m  long  optical  fiber 
so  as  to  modulate  the  uncorrelated  background  noise  in  the  tail  of  the  correlation  function. 
Here,  we  measure  ffcohd(0)  =  1-2.  In  the  second  experiment  we  apply  modulation  at  random 
times,  using  an  external  10  MHz  digital  signal  as  a  trigger  for  the  function  generator.  As 
expected,  we  observe  a  reduced  rate  of  paired  counts  and  no  change  in  the  shape  of  the 
correlation  function. 

The  method  demonstrated  in  this  work  might  be  used  to  optimally  load  a  single  photon 
into  an  optical  cavity,  or  instead,  to  study  the  transient  response  of  atoms  to  different  single 
photon  waveforms.  In  the  context  of  light-matter  interfaces,  it  may  improve  the  efficiency 
of  storage  and  retrieval  of  single  photons  in  atomic  ensembles.  For  quantum  information 
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applications,  both  amplitude  and  phase  modulators  could  be  used  to  allow  full  control  over 
the  single  photon  waveforms.  For  example,  one  could  a  construct  a  single  photon  waveform 
that  is  a  train  of  identical  pulses  with  information  encoded  into  the  relative  phase  differ¬ 
ence  between  consecutive  pulses.  The  importance  of  the  electro-optic  method  is  its  speed 
and  ability  to  modulate  phase  as  well  as  amplitude.  The  technique  provides  the  technol¬ 
ogy  for  studying  the  response  of  atoms  to  shaped  single  photon  waveforms  on  a  time-scale 
comparable  to  the  natural  linewidth. 

For  more  details  please  see:  P.  Kolchin,  Chinmay  Belthangady,  Shengwang  Du,  G.Y.  Yin, 
and  S.E.  Harris,  “Electro-Optic  Modulation  of  Single  Photons,”  Phys,  Rev.  Lett.  101, 
103601  (September,  2008). 


5  Modulation  and  Measurement  of  Time-Energy  En¬ 
tangled  Photons 

In  2008  we  reported  the  demonstration  of  conditional  temporal  shaping  of  single  photon  wave¬ 
forms  using  electro-optic  modulators.  In  this  report,  we  extend  this  work  to  the  modulation 
of  biphotons.  A  highlight  is  the  invention  and  first  experimental  demonstration  of  a  Fourier 
technique  that  allows  the  measurement  of  fast  biphotons  using  slow  detectors.  Figure  1  shows 
(a)  the  schematic  of  our  proof-of-principal  experiment  and  (b)  the  laser /rubidium- atom  in¬ 
teraction  diagram  for  paired  phoion  generation.  Parametrically  down-converted  spontaneous 
signal  and  idler  photons,  or  as  in  the  experiment  of  this  work,  Stokes  and  anti-Stokes  pho¬ 
tons,  are  incident  on  synchi*onoiisly  driven  sinusoidal  amplitude  modulators.  Without  any 
modulation,  in  general,  the  setup  will  directly  yield  the  Glauber  correlation  function 
where  r  is  the  relative  arrival  time  of  the  signal  and  idler  photons,  with  its  time  resolution 
limited  by  the  speed  of  photon  detector  (SPCM)  and  TDC  system.  As  shown  below,  our 
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new  method,  by  adding  AM  photon  modulation,  will  give  temporal  biphoton  correlation 
measurement  with  time  resolution  beyond  the  SPCM-TDC  limit. 

When  amplitude  modulators  mi{t)  and  m2{t)  are  introduced  between  the  down  conversion 
source  and  the  detectors,  the  modulated  correlation  function  is  then  written  as 


(5.1) 


where  the  subscripts  Mand  0  indicate  cases  of  modulation  on  or  off.  In  real  experiments, 
photons  arrive  at  a  random  time  t  which  is  averaged  out  over  the  data  collection  period  T: 


=  \Tni{t)f  \Tn2{t  +  T)f  dt, 


(5.2) 


where  M.{t)  is  the  intensity  (Correlation  function  of  the  modulators  in  the  signal  and 
idler  channels.  If  both  channels  are  modulated  by  sinusoidal  amplitude  modulators  with 
frequency  u  and  a  common  phase  (p,  that  is  by  modulators  m\{t)  =  m2{t)  =  cos(a;t  +  (p), 
then,  irrespective  of  this  phase,  M[t)  =  1/4  +  l/8cos(2a;r). 

If  the  detectors  are  slow  in  the  sense  that  they  integrate  over  the  length  of  the  bipho¬ 
ton  wave  packet  ,  but  short  as  ( ompared  to  the  inverse  rate  of  bi photon  generation,  the 
measurement  becomes  an  integral 

poc - poo 

/  G^^(r,C(;)dr  =  1/8  /  [2  +  cos(2c(;t)]  Go^^(T)(iT.  (5.3) 

Jo  Jo 

We  neglect  the  DC  term  and  normalize  to  obtain  the  Fourier  cosine  transform  pair 

F{2u’)  =  •/—  J  G^q\t)  cos{2L0T)dT, 

Pi 

=  F{2L0)cos{2cuT)duj.  (5.4) 

In  the  measurement  procedure,  F{2uj)  is  the  measured  coincidence  count  rate  between  low 
speed  single  photon  counting  modules  (SPCMs)  as  a  function  of  the  sinusoidal  modulation 
frequency  ui.  The  slow  detection  system  (SPCM+TDC)  means  that  hardware  integrates 
signals  over  r,  which  is  the  relative  arrival  time  of  the  signal  and  idler  photons.  We  then 
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Figure  7:  (a)  Schematic  of  experiment.  Long  biphotons  are  generated  in  the  cold-atom  Rb 
cell.  Using  an  optical  fiber,  the  anti-Stokes  photon  is  delayed  by  175  ns,  and  the  photons  are 
modulated  by  synchronously-driven  sinusoidal  modulators  before  correlation  detection  with 
a  time-to-digital  converter,  (b)  Energy  level  diagram  for  paired  photon  generation  in  Rb. 

apply  the  above  inverse  Fourier  cosine  transform  to  yield  the  Glauber  correlation  function 
G^\t)  and  therefore  the  square  of  the  absolute  value  of  the  biphoton  wavefunction. 

Our  experiment  makes  use  of  long  biphotons  that  are  produced  using  the  techniques  of 
electromagnetically  induced  transparency  and  slow  light.  The  use  of  long  biphotons  allows 
us  to  compare  the  correlation  function  measured  by  our  Fourier  transform  technique  with 
a  direct  measurement  using  fast  detectors  and  a  TDC.  The  experimental  configuration  is 
shown  in  Fig.  7(a).  Paired  photons  are  generated  with  cold  Rb  atoms  using  the  method 
of  Balic  et  al.  We  apply  strong  counterpropagating  pump  and  coupling  lasers  (not  shown) 
to  produce  phase  matched  counter-propagating  pairs  of  time-energy  entangled  Stokes  and 
anti-Stokes  photons.  We  use  a  ^''^Rb  two-dimensional  magneto-optic  trap  with  an  optical 
depth  which  can  be  varied  between  10  to  60  to  generate  biphotons  with  temporal  lengths 
between  50  and  900  ns.  The  inset  in  Fig.  8(a)  shows  the  biphoton  wavefunction  obtained  at 
an  optical  depth  of  35  as  directly  measured  using  a  TDC.  Two  features  are  of  interest.  First, 
the  width  of  the  wavefunction  is  determined  by  the  slow  group  velocity  of  the  anti-Stokes 
photon  and  varies  linearly  with  the  optical  depth.  Second,  the  distinctive  sharp  feature  at 
the  leading  edge  is  a  Sommerfield-Brillouin  precursor  that  ensures  that  the  earliest  signal 
reaches  a  detector  at  the  speed  of  light  in  vacuum. 
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Figure  8:  (color  online)  Modulation  of  the  biphoton  correlation  function,  (a)  mi{t)  =  1  and 
m2{t)  =  cos{cjt  +  if).  The  inset  in  (a)  shows  the  correlation  function  with  both  modulators 
open,  (b)  7ni{t)  =  m2(t)  =  cos(ut  4“  ^).  Here  u  =  27rx  35  xlO^. 

The  generated  Stokes  and  anti-Stokes  photons  are  transmitted  through  10  GHz  electro¬ 
optic  amplitude  modulators  (Eospace  Inc.)  with  a  half-wave  voltage,  14  of  1.3V.  To  obtain 
a  perfectly  sinusoidal  output,  the  modulators  are  biased  at  maximum  transmission  and  the 
input  voltage  is  varied  linearly  using  a  triangular  waveform  that  varies  between  —14  and  +14. 
This  waveform  is  generated  by  a  fast  function  generator  (Tektronics  AFG3252)  with  two 
output  channels  whose  frequencies  and  phases  can  be  varied  independently.  The  modulated 
photons  are  then  sent  to  SPCMs  (Perkin  Elmer  SPCM-AQR-13),  which  are  connected  to  the 
start  and  stop  inputs  of  a  TDC  (Fast-Comtec  TDC  7886S).  Coincidence  counts  are  binned 
into  histograms  and  plotted  as  a  function  of  the  time  difference  between  the  detection  of  a 
Stokes  and  an  anti-Stokes  photon. 

In  Fig.  8  (color  online)  we  demonstrate  the  modulation  of  a  biphoton  wavefunction.  The 
data  are  recoded  by  binning  coincidence  counts  versus  time  into  1  ns  bins.  In  Fig.  8(a),  the 
Stokes  modulator  is  turned  off  and  is  biased  at  maximum  transmission;  and  the  anti-Stokes 
modulator  is  driven  at  35  MHz.  In  agreement  with  theory,  the  biphoton  wavefunction  is  not 
modulated.  In  Fig.  8(b),  both  modulators  are  modulated  at  35  MHz  with  the  same,  but  ar¬ 
bitrary,  phase.  The  correlation  function  is  now  modulated  at  twice  the  applied  frequency,  i.e. 
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at  70  MHz.  When  driven  by  non-siriusoidal  waveforms,  the  correlation  function  is  modulated 
by  the  cross-correlation  of  the  two  modulating  signals.  In  agreement  with  Eq.  (5.2),  we  have 
verified  that  when  the  two  modulating  signals  are  square  waves  with  the  same  frequency,  the 
correlation  function  is  modulated  by  a  triangular  function. 


Figure  9:  Fourier  Transform  measurement  technique  at  an  optical  depth  of  15(Left)  and 
35(Right).  Upper:  Frequency  domain  data.  Lower:  Fourier  Cosine  transforms  (red)  of 
corresponding  frequency  domain  data,  and  real-time  fast  direct  temporal  correlation  data 
(blue).  The  term  modulation  frequency  on  the  x-aixis  of  (a)  and  (c)  refers  to  the  applied 
frequency.  The  observed  modulation  frequency  [Eq.  (5.4)]  is  a  factor  of  2  higher. 

We  next  demonstrate  the  Fourier  measurement  technique.  A  35  m  long  polarization- 
maintaining  fiber  is  used  to  delay  the  anti-Stokes  photon  by  175  ns  in  the  experiment.  The 
modulators  are  driven  synchronously  at  frequencies  between  0  and  30  MHz.  To  simulate 
slow  detectors,  we  bin  coincidences  into  1  /is  bins,  and  since  the  temporal  extent  of  the 
waveform  is  less  than  1  /iS,  the  entire  biphoton  is  contained  within  the  first  time  bin.  From 
the  Fourier  transform  property  that  delay  in  the  time  domain  corresponds  to  oscillation  in 
the  frequency  domain,  the  plot  of  coincidence  counts  versus  modulation  frequency  shows 
ripples  with  a  frequency  (5.7  MHz)  equal  to  the  inverse  of  the  time  delay.  We  show  two 
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data  sets  in  Fig.  9  for  optical  depth  of  15  (left)  and  35  (right),  for  relatively  short  and 
long  biphoton  wave  packets,  respectively.  In  Figs.  9(b)  and  (d),  we  compare  the  biphoton 
wavefunction  as  directly  measured  with  1  ns  bins  (blue  line)  to  the  Fourier  cosine  transform 
of  the  traces  in  Figs.  9(a)  and  (c),  (red  line),  vertically  scaled  to  match  the  peak  value.  We 
find  reasonable  agreement  between  the  two  methods.  We  note  that  the  sharp  spike  in  the 
temporal  trace  near  t  =  0  results  from  the  DC  component  in  the  frequency  domain  trace. 

We  have  shown  how  biphotons  may  be  modulated,  and  how  this  modulation  may  be  used 
to  measure  the  magnitude  of  the  biphoton  wavefunction.  Though  we  have  used  long  bipho¬ 
tons  and  low  modulation  frequenci(3s,  we  believe  that  this  technique  should  be  extendable 
to  short  biphotons.  A  commercially  available  telecommunication  modulator  driven  at  a  fre¬ 
quency  of  60  GHz  and  therefore  modulating  at  120  GHz  will  allow  measurement  of  biphotons 
with  a  minimum  length  of  about  8  picoseconds.  This  is  about  a  factor  of  five  faster  than 
state-of-the-art  commercial  SPCMs.  If  the  comparison  is  made  on  the  basis  of  state-of-the- 
art  polymer  light  modulators  operating  at  200  GHz  and  therefore  modulating  at  400  GHz, 
then  the  Fourier  technique,  will  allow  measurement  down  to  about  2.5  ps.  This  is  about 
eight  times  faster  than  the  fastest  reported  superconducting  detector.  Looking  further  to 
the  future,  all-optical  light  modulators  have  been  demonstrated  at  frequencies  greater  than 
1  THz;  thereby  in  principle  allowing  measurement  of  photons  on  femtosecond  time  scales. 
We  note  that  an  advantage  of  the  Fourier  technique  as  compared  to  either  sum  frequency 
correlation  or  Hong-Ou-Mandel  interference  is  that  the  signal  and  idler  photons  may  be  cor¬ 
related  at  distant  detectors  and  do  not  need  to  be  brought  together  at  a  summing  crystal  or 
at  a  beam  splitter. 

For  more  details  please  see:  Chinmay  Belthangady,  Shengwang  Du,  Chih-Sung  Chuu,  G.Y. 
Yin,  and  S.E.  Harris,  “Modulation  and  measurement  of  time-energy  entangled  photons”, 
Phys.  Rev.  A,  Rapid  Communications,  80,  031803  (September  2009). 


21 


6  Resonator  Sum  Frequency  Generation  with  Time- 
Energy  Entangled  Photons 

We  recently  developed  a  technique  for  enhancing  the  quantum  component  of  sum  frequency 
radiation  that  is  obtained  by  summing  pairs  of  signal  and  idler  photons.  By  resonating  the 
sum  frequency  we  observe  generation  that  varies  linearly  with  input  power  and  is  increased 
by  a  factor  of  twelve.  The  essence  of  the  technique  is  similar  to  the  Purcell  effect  where  the 
work  done  by  a  generated  dipole  moment,  in  this  case  monochromatic,  against  an  electric 
field  at  its  own  frequency  is  increased  by  the  presence  of  the  resonator.  Of  importance, 
this  resonance  enhancement  occurs  pair  by  pair,  thus  resulting  in  the  linear  dependence  of 
generated  sum  power  on  the  incident  biphoton  rate. 

The  technique  of  resonant  enhancement  is  shown  in  Fig.  10.  Time-energy  entangled 
photons  termed  as  the  signal  and  idler  are  generated  by  a  parametric  down-converter  (not 
shown)  with  a  monochromatic  pump  so  that  the  pump  frequency  ujp  equals  the  sum  of  the 
signal  and  idler  frequencies.  These  photons  are  summed  in  a  nonlinear  crystal  placed  inside 
a  one-sided  cavity  whose  longitudinal  mode  at  frequency  ujq  is  detuned  from  the  pump 
frequency  by  5u)q  =  ujp  —  uOq. 

Rl=1  R2=R 


crystal 

Figure  10:  Resonant  enhancement  technique. 

The  SFG  rate  consists  of  two  terms,  a  term  due  to  quantum  SFG,  and  a  term  due  to 
classical  SFG.  It  can  be  shown  that  the  quantum  term  is  monochromatic  and  varies  linearly 
with  input  power,  while  the  classical  term  has  a  lineshape  equal  to  the  convolution  of  the 
signal  and  idler  spectra  and  varies  quadratically  with  input  power.  Similarly,  the  SFG  rate  of 
the  traveling  wave  case  (without  cavity)  also  consists  of  a  quantum  term  that  varies  linearly 
with  input  power  and  a  classical  term  that  varies  quadratically. 
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The  quantum  SFG  rate  is  enhanced  by  the  presence  of  the  cavity  by  an  amount  equal  to 
the  ratio  of  the  quantum  terms  for  the  two  cases.  When  the  cavity  is  tuned  to  resonate  the 
pump  frequency  {ScOq  =  0),  the  enhancement  is 

.  /  4(1 -i?) 

cavity /traveling  = 

The  largest  quantum  SFG  enhancement  is  obtained  for  i?  =  1  —  2^;  for  this  reflectivity  the 
enhancement  ratio  is  equal  to  the  inverse  of  the  round  trip  power  loss,  l/(2<^). 

Of  importance,  assuming  that  the  signal  and  idler  bandwidths  are  large  as  compared  to 
the  cavity  free  spectral  range  A,  enhancement  by  the  use  of  a  resonant  cavity  is  only  possible 
in  the  quantum  and  not  in  the  classical  regime.  By  using  the  resonance  technique  described 
here,  both  the  quantum  SFG  rate  and  the  ratio  of  the  quantum  to  classical  term  may  be 
improved.  Our  proof-of-principal-experiment  configuration  is  shown  in  Fig.  11.  Frequency- 


Figure  11:  Schematic  of  experiment. 

degenerate  biphotons,  with  a  calculated  bandwidth  of  32  nm,  are  generated  by  spontaneous 
parametric  down-conversion  in  a  20  mm  long  periodically-poled,  magnesium  oxide-doped, 
stoichiometric  lithium  tantalate  crystal  (MgO:PPSLT)  pumped  by  an  8  W,  CW  laser  at  532 
nm.  A  four-prism  setup  is  used  to  filter  out  the  strong  pump  and  to  provide  dispersion 
compensation  for  the  down-converted  biphotons.  SFG  of  these  biphotons  occurs  inside  a 
cavity  that  is  reflecting  at  532  nm  and  transparent  at  1064  nm. 

The  SFG  crystal  is  a  20  mm  long  MgOiPPSLT  crystal  anti-reflection  coated  at  1064  nm 
and  532  nm.  The  confocal  cavity  consists  of  two  20  mm  radius  of  curvature  spherical  mirrors, 
one  of  which  is  mounted  onto  a  piezoelectric  crystal.  The  input  mirror  has  a  reflectance  of 
Rin  >  99.5%,  and  the  output  mirror  has  a  reflectance  of  Rout  =  95.5%.  (The  measured 
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single-pass  crystal  loss  is  ^  ~  2%,  and  the  optimum  reflectivity  of  the  output  mirror  is 
=  1  —  2^  ~  96%).  The  SFG  signal  is  filtered,  then  collected  into  multimode  fiber  and 
detected  with  a  single  photon  counting  module. 

A  typical  scan  of  the  cavity  over  a  free  spectral  range  is  shown  in  Fig.  12(A).  The  resonant 
behavior  of  the  SFG  process  is  clearly  observed,  with  peak  generation  rates  on  resonance 
that  are  significantly  higher  than  the  traveling  wave  rate  of  3200  counts/s.  While  the  cavity 
is  passively  stabilized  to  within  a  fraction  of  a  free  spectral  range  during  the  scan  period, 
temperature  instability  and  air  currents  cause  the  peak  measured  SFG  rate  to  fluctuate  by 
an  amount  that  is  6  times  larger  than  the  shot  noise. 

To  better  quantify  the  enhancement,  we  measure  the  resonantly-enhanced  SFG  rate  at 
various  down-converted  infrared  (IR)  input  powers.  Figure  12(B)  shows  the  averages  and 
standard  deviations  of  the  peak  SFG  rates  at  each  input  power.  Also  shown  is  the  traveling- 
wave  SFG  rate  obtained  by  removing  the  out-coupling  mirror.  The  linear  fits  to  the  data  in 
Fig.  12(B)  demonstrate  two  results.  First,  the  ratio  of  their  slopes  gives  the  SFG  enhance¬ 
ment  ratio  of  12.  Second,  the  linear  dependence  of  the  SFG  rate  on  input  power,  in  both 
cases,  shows  that  it  is  the  quantum  SFG  term  that  is  enhanced.  The  enhancement  is  a  factor 
of  two  less  than  the  theoretical  prediction.  The  discrepancy  may  be  due  to  imperfect  mode 
matching  or  cavity  instability. 

In  addition  to  demonstrating  the  effect  experimentally,  we  have  also  developed  a  theory 
for  resonantly-enhanced  quantum  SFG.  We  work  in  the  Heisenberg  picture  with  traveling 
wave  signal  and  idler  beams  denoted  by  the  annihilation  operators  as{t,  z)  = 
ds{t^  z)  exp  [—i{uJst  —  kgz)]  and  ai(t,  z)  =  di{t^  z)  exp  [—i{uJit  —  kiz)]^  where  operators  with  a 
tilde  vary  slowly  with  time  and  distance.  We  use  a  slowly  varying  envelope  formalism  with 
the  standing-wave  cavity  mode  operator  written  as  aq{t^  z)  =  dq{t)  exp{—iujqt)  sm{kqz)^  where 
kq  =  qn/L.  We  project  the  broadband  generated  dipole  moment  operator,  proportional  to 
ds{t,z)di{t^z)^  against  the  cavity  mode.  With  the  traveling  wave  SFG  field  emitted  from 
mode  q  denoted  by  Usum(^j^)  =  a5um(^)  ^xp  [— i(cjpt  —  fcpz)],  the  equation  for  the  evolution 
of  dq{t)  and  its  relation  to  the  SFG  field  is: 
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Figure  12:  (A)  The  cavity  is  scanned  over  a  free  spectral  range,  thereby  demonstrating 
the  SFG  resonance  effect.  The  pink  solid  line  is  the  theoretical  plot  of  the  quantum  term 
normalized  to  the  average  peak  rate.  Blue  circles  are  experimental  data.  (B)  SFG  rate  for 
cavity  and  traveling  wave  experiments.  Circles  are  resonant  SFG  rates  averaged  over  20-25 
peaks  obtained  in  slow  cavity  scans.  Squares  are  traveling-wave  SFG  rates  averaged  over 
60-second  intervals. 
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(6.1) 


+  ^aq{t)  =  Pq{t)  exp(-z(5wi)  +  P{t) 


Denoting  the  spacing  of  the  cavity  modes  at  the  sum  frequency  by  A  =  c/ (2Ln),  the  decay 


rate  F  of  photons  in  the  cavity  is  determined  by  the  mirror  reflectivity  R  and  the  single  pass 
power  loss  With  the  output  coupling  rate  7  =  A(1  —  R),  the  total  (power)  decay  rate  is 
r  =  2^A  +  7.  We  take  all  fields  to  be  plane  waves  with  cross  sectional  area  A  and  take  the 
refractive  index  n  and  nonlinearity  d  to  be  independent  of  frequency.  With  the  cavity  and 
summing  crystal  of  the  same  length,  the  coupling  constant  =  {d/n'^){fXohLOpLOstOiL/ . 
The  quantity  P{t)  is  a  Langevin  noise  operator  that  has  contributions  from  an  incoming 
wave  at  the  right  hand  mirror  (not  shown)  and  a  macroscopic  loss  term.  Both  are  negligible 
at  room  temperature.  The  normalization  is  such  that  {a)q{t)dq{t))  is  the  total  number  of 
photons  in  the  cavity  mode,  and  {al^j^(t)asum(t))  is  the  rate  of  sum  photons  exiting  the 


cavity. 

The  solution  of  Eqs.  (1)  for  the  SFG  field  operator  is 


(6.2) 


We  assume  that  the  signal  and  idler  fields  are  not  depleted  in  the  summing  crystal  and 
transform  to  the  frequency  domain  using  as{t,z)  =  f  as(LJ,  z)  exp(—iu;t)doj  and  ai(t,z)  = 


f  ai{uj,z)exp{—iujt)(Lj.  Pq{t)  becomes 


X  exp  +U)2  —  ojp)t\  du!iduJ2 


^{u)i,u)2)  =  exp  ~i 


=  exp  -~i - ^ -  sine 


Afc(a;i,ci;2)Tl  .  rAfc(a;i,a;2)T’ 


2  2 


(6.3) 
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where  $(a;i,6(;2)  is  a  phase  matching  factor  with  sinc(a;)  =  sin(x)/x  and  Ak{u}i,L02)  =  kq  — 
[A:(t<;i)  +  A;(a;2)]-  Equation  (6.3)  is  substituted  into  Eq.  (6.2)  to  obtain  the  SFG  output  field: 


asum[  )  -  J_^  -  iiuj,  +u;2-u;,  +  6cv) 

X  exp  [—i{ijO\  +  LU2  —  <^p  +  du\dw2-  (6.4) 

For  monochromatically-pumped  down  conversion,  the  down-converted  signal  and  idler 
fields  at  the  input  of  the  summing  crystal,  0^(0;)  and  aj(a;'),  where  to'  =  Up  —  oj,  can  be 
described  in  terms  of  initial  vacuum  fields  aso(u^)  and  0^0(0;')  as 

as(u;)  =  A(u;)aso(c^)  +  5(a;)a|'o(a;') 

a|(a;')  =  C(a;)a5o(a;)  +  F)(a;)ato(a;').  (6.5) 

Substituting  Eqs.  (6.5)  into  Eq.  (6.4)  and  noting  the  commutator  [aio(wi),  aJo(‘*^2)]  = 
(l/27r)^ij^(u;i  —  a;2),  we  evaluate  (cisum(^)^sum(i))  and  find  the  rate  of  SFG  photons  in  mode 
q  exiting  the  cavity  to  be 


T^r 


-  (?)7: 


7 


F2  +  ASlu^ 


f{uj)  + 


7 


F2  +  4(a;  -  a;,)2 


duj, 


where  f{uj)  and  ^(a;)  are 


fi<^) 

9(1^) 


/oo 

A{n)c*{n)^n,u;p-n)dn 

■00 

/oo 

\B{n)f\C{n-u;  +  u;p)fmn 

■00 


5(a;  —  ojp) 

,u;-n)fdn. 


(6.6) 


(6.7) 


The  SFG  rate  in  Eq.  (6.6)  consists  of  a  quantum  term  containing  the  function  /(a;)  and  a 
classical  term  containing  the  function  g{u;)  .  The  quantum  term  is  the  result  of  generation 
with  correlated  photons;  it  may  be  shown  to  scale  linearly  with  input  power  and  has  a 
monochomatic  spectrum.  The  classical  term  comes  from  SFG  with  uncorrelated  photons 
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and  varies  qiiadratically  with  input  power.  Its  spectrum  is  proportional  to  the  convolution 
of  the  signal  and  idler  spectra,  multiplied  by  the  cavity  lineshape. 

In  order  to  normalize  the  improvement  that  results  from  the  use  of  a  cavity,  we  write  the 
SFG  rate  for  the  traveling  wave  case  with  no  cavity  present.  This  rate  is 


/  .  2 

T^tw  =  j  J  [/(^)  +  9{(^)]  du;,  (6.8) 

where  Ktw  =  and  kq  in  the  function  $(cc;i,a;2)  becomes  k{u;\  +002)-  As  in  the  cavity 

case,  the  SFG  rate  consists  of  a  quantum  term  that  varies  linearly  with  input  power  and  a 
classical  term  that  \^ries  quadratically. 


The  quantum  SFG  rate  is  enhanced  by  the  presence  of  the  cavity  by  an  amount  equal  to 
the  ratio  of  the  terms  containing  /(a;)  in  Eqs.  (6.6)  and  (6.8).  When  the  cavity  is  tuned  to 
resonate  the  pump  frequency  (5c<;  =  0),  the  quantum  SFG  enhancement  ratio  is 


^  4(1 -i?) 

"  [(l-i?)  +  2^f 


(6.9) 


Maximizing  with  respect  to  i?,  the  largest  quantum  SFG  enhancement  is  obtained  for  R  = 
1  —  2G  for  this  reflectivity  the  enhancement  ratio  is  equal  to  the  inverse  of  the  round  trip 
power  loss, 


To  study  the  effect  of  the  cavity  on  the  classical  component  of  SFG,  we  assume  that 
the  signal  and  idler  spectra  are  sufficiently  broadband  that  the  dipole  moment  at  the  sum 
frequency  is  constant  over  many  cavity  modes.  In  this  case,  g{ijo)  is  a  constant,  and  the 
generated  SFG  spectrum  is  periodic  with  resonant  peaks  separated  from  each  other  by  the 
cavity  free  spectral  range.  To  obtain  the  net  enhancement  of  the  classical  term,  we  integrate 
over  a  single  free  spectral  range  centered  at  Uq.  The  resulting  classical  SFG  enhancement 
ratio  is  then 
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(6.10) 
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where  the  second  equality  follows  for  sufficiently  high  finesse  that  Tan“^(27rA/r)  — >  7r/2. 
If  the  cavity  is  lossless,  r/c  =  1.  Though  generation  at  resonance  is  enhanced,  generation 
off-resonance  is  suppressed  so  that  the  integrated  classical  SFG  enhancement  ratio  is  unity. 

Equations  (6.9)  and  (6.10)  show  how  resonating  the  SFG  field  enhances  the  efficiency  of 
quantum  SFG  while  leaving  classical  SFG  unchanged.  The  enhancement  of  quantum  SFG 
may  find  use  in  correlation  experiments  that  measure  the  lengths  of  biphoton  wavefunctions 
which  are  too  short  to  be  measured  with  time-resolved  detectors.  In  such  experiments,  the 
SFG  crystal  must  be  very  short  to  have  a  sufficiently  large  acceptance  bandwidth  for  the 
biphotons,  and  the  SFG  signal  may  be  small  as  a  result.  A  resonator  could  enhance  the 
signal  to  measurable  levels. 


7  Observation  of  Nonlocal  Modulation  with  Entangled 
Photons 

The  idea  of  what  we  term  here  as  nonlocal  modulation  was  published  about  one  year  ago; 
i.e.,  S.  E.  Harris,  “Nonlocal  Modulation  of  Entangled  Photons,”  Phys.  Rev.  A  78,  021807(R) 
(2008).  This  section  of  the  report  describes  the  first  experimental  observation  of  this  effect. 

When  the  photons  of  a  time-energy  entangled  pair  are  sent  through  different  channels 
having  arbitrary  dispersions,  the  dispersion  in  one  channel  may  be  negated  by  dispersion  of 
the  opposite  sign  in  the  other  channel.  This  effect  results  from  a  quantum  mechanical  inter¬ 
ference  and  has  no  classical  analog.  This  is  now  termed  as  nonlocal  dispersion  compensation. 
We  report  the  first  observation  of  a  time-frequency  analog  to  nonlocal  dispersion  cancella¬ 
tion  and  term  this  effect  as  nonlocal  modulation.  Consider  a  simplified  concept  system  as 
shown  in  left  portion  of  Fig.  13,  where  a  monochromatic  pump  generates  non-degenerate 
time-energy  entangled  photon  pairs.  The  signal  and  idler  photons  (Channels  1  and  2)  pass 
through  sinusoidal  phase  modulators.  These  modulators  are  driven  at  the  same  modulation 
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frequency,  and  their  relative  phase  may  be  varied.  After  passing  through  the  modulators, 
the  signal  and  idler  photons  are  dispersed,  for  example,  by  a  prism,  and  the  relative  po¬ 
sitions  (frequencies)  of  the  signal  and  idler  photons  are  correlated.  When  the  modulation 
frequency  is  small  as  compared  to  the  spectral  bandwidth  of  the  signal  or  idler,  we  find  a 
consequence  of  time-energy  entanglement  that  we  term  as  nonlocal  modulation.  Specifically, 
these  distant  modulators  act  cumulatively  to  determine  the  apparent  modulation  depth.  If 
the  two  identical  modulators  have  opposite  phase,  they  negate  each  other  and  act  as  if  nei¬ 
ther  modulator  were  present.  Conversely,  if  operated  with  the  same  phase,  they  produce  the 
same  correlation  as  does  a  single  modulator  with  twice  the  modulation  depth  acting  on  only 
one  of  the  photons. 


Figure  13:  Nonlocal  modulation.  Signal  and  idler  photons  are  phase  modulated  at  the  same 
frequency  and  with  controllable  phase.  The  signal  and  idler  beams  are  frequency-dispersed, 
and  the  positions  (frequencies)  of  the  detected  photons  on  the  photodetectors  are  correlated. 
Left:  simplified  concept  diagram.  Right:  experimental  configuration. 

If  the  resolution  of  the  frequency  correlator  of  Fig.  13  (left)  is  infinite,  the  resulting 
correlation  traces  for  phased  and  anti-phased  modulators  may  look  something  like  Fig.  14. 
Since  the  pump  laser  which  generates  the  entangled  photons  is  monochromatic,  they  are 
delta-function  correlated  in  frequency  before  any  modulation  occurs.  The  modulators  act 
on  this  delta  function  correlation  by  generating  sidebands  whose  distribution  depends  on 
the  amplitude  of  the  sinusoidal  driving  signal.  If  the  modulators  are  run  in  phase  with 
each  other,  the  distribution  of  sidebands  resembles  that  of  a  single  modulator  with  twice 
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Figure  14:  Example  frequency  correlation  function  with  (a)  phased  and  (b)  anti-phased 
sinusoidal  modulators.  Because  the  modulators  have  the  same  amplitude,  cumulative  mod¬ 
ulation  results  in  either  (a)  doubling  of  modulation  depth  or  (b)  cancellation  of  modulation 
in  the  spectral  correlation. 

the  driving  amplitude  [Fig.  14(a)].  If  run  with  opposite  phase,  the  modulators  cancel  each 
other,  and  a  single  delta  function  is  present  in  the  correlation,  as  if  neither  modulator  were 
present  [Fig.  14(b)]. 

The  more  detailed  experimental  configuration  is  shown  in  right  part  of  Fig.  13.  We 
pump  a  20  mm  long,  periodically-poled,  magnesium  oxide-doped  stoichiometric  lithium  tan- 
talate  crystal  (PPSLT)  with  0.8  W  from  a  532  nm  cw  laser.  The  nonlinear  crystal  is  phase 
matched  to  produce  32  nm  bandwidth,  degenerate  photon  pairs  at  1064  nm.  All  fields  are 
polarized  along  the  extraordinary  axis  of  the  crystal.  The  generated  photons  are  filtered 
from  the  strong  532  nm  pumping  beam  using  a  four-prism  setup  and  are  then  coupled  into  a 
polarization-maintaining  fused-fiber  beam  splitter  which  diverts  the  photons  into  Channels 
1  and  2  with  equal  probability.  The  photons  are  passed  through  identical  sinusoidal  phase 
modulators  (EOSPACE)  which  are  driven  at  30  GHz  with  modulation  depths  of  about  1.5 
radians.  The  relative  phase  between  the  modulators  is  controlled  using  a  calibrated  phase 
trimmer.  The  photons  then  pass  through  identical  monochromators,  each  having  a  linear 
dispersion  of  210  GHz/mm  and  a  Gaussian  instrument  response  function  with  a  FWHM  of 
8.5  GHz.  To  obtain  frequency  domain  correlation  curves,  we  fix  the  output  slit  in  Channel 
1  at  Xi  and  scan  the  position  X2  in  Channel  2.  The  photons  that  are  transmitted  through 
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the  monochromator  slits  are  coupled  into  multimode  fibers  and  detected  with  time-resolved 
single  photon  counting  modules  (SPCMs,  id  Quantique  id400  and  PerkinElmer  SPCM-AQR- 
16-FC). 

The  primary  experimental  results  of  this  work  are  shown  in  Fig.  15.  For  each  case,  we  set 
the  monochromator  slit  in  Channel  1  at  an  arbitrary  position  xi  that  is  near  the  center  of 
the  generated  32  nm  spectrum  and  leave  the  position  of  this  slit  fixed  thereafter.  The  slit  in 
Channel  2  is  scanned  over  positions  X2^  and  the  coincidence  rate  of  the  two  detectors  (with 
gate  width  T  =  1.25  ns)  is  recorded  as  a  function  of  this  position.  For  each  position,  the 
rate  is  averaged  for  20  s. 

With  the  pump  frequency  defined  as  and  the  position  X2  proportional  to  the  frequency 
iJ2^  we  express  the  coincidence  rate  as  a  function  of  relative  frequency  A  =  cj2  —  {cjp  —  Ui). 
The  scale  of  the  frequency  axis  is  calibrated  by  measuring  the  sideband  spacing  of  a  single¬ 
mode  1064  nm  laser  modulated  at  30  GHz,  with  the  zero  position  chosen  (at  the  start  of  the 
experiment)  as  the  location  of  the  correlation  peak  for  unmodulated  photon  pairs. 

The  upper  portion  of  Fig.  15  shows  the  experimental  results  without  modulation  and 
with  modulation  in  a  single  channel.  In  Fig.  15(a)  (upper),  both  modulators  are  turned 
off  by  disconnecting  their  30  GHz  drive  signals.  As  expected  by  energy  conservation,  a 
single  correlation  peak  is  observed.  In  Fig.  15(b)  (upper),  Channel  1  is  phause  modulated 
as  exp[i(5sin(a;^t)]  with  a  modulation  depth  of  5  =  1.5,  and  Channel  2  is  not  modulated. 
The  frequency  correlation  is  now  distributed  over  a  set  of  sidebands,  having  Bessel  function 
amplitudes  J^(^),  whose  total  area  is  equal  to  that  of  upper  (a)  of  Fig.  15. 

In  Fig.  15(a)  (lower),  both  modulators  are  turned  on  at  a  modulation  depth  ofS  =  1.5,  and 
the  cable  length  is  adjusted  so  that  they  have  the  same  phase.  They  now  act  cumulatively 
(constructively  interfere)  to  produce  a  set  of  sidebands  having  a  Bessel  function  distribution 
J^{2S).  The  frequency-domain  correlation  function  of  two  distant  modulators  is  therefore 
the  same  as  that  which  would  be  obtained  by  correlating  an  unmodulated  photon  with  a 
photon  modulated  at  twice  the  modulation  depth. 

In  Fig.  15(b)  (lower),  the  modulators  are  run  at  the  same  depth  as  in  the  previous  para¬ 
graph,  but  now  the  relative  cable  length  is  adjusted  so  that  the  modulators  are  run  in  phase 
opposition.  The  modulators  now  destructively  interfere,  and  no  sidebands  are  visible. 
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Figure  15:  Frequency  correlation  measurements.  The  horizontal  axes  measure  relative  fre¬ 
quency  changes  as  slit  position  X2  is  scanned.  Upper:  (a)  both  modulators  turned  off  and 
(b)  the  modulator  in  Channel  1  running  at  a  modulation  depth  of  1.5.  Lower:  both  modu¬ 
lators  running  (a)  with  the  same  phase  and  (b)  with  opposite  phase.  The  modulation  depth 
for  both  modulators  is  1.5.  Circles  are  data;  curves  are  theoretical  fits  (see  text). 


The  striking  difference  between  Figs.  14  and  15  is  that  the  sidebands  are  not  delta  func¬ 
tions  but  have  widths  determined  by  the  resolution  of  the  frequency  correlator  in  Fig.  13 
(right).  We  have  developed  a  theory  to  predict  the  measured  correlation  functions  with 
finite-resolution  monochromators  as  used  in  our  experiment.  The  following  describes  this 
theory,  which  is  used  to  produce  the  solid  curves  in  Fig.  15. 

Working  in  the  Heisenberg  picture,  a  nonlinear  crystal  of  length  L  is  pumped  by  a 
monochromatic  laser  at  frequency  Up.  A  positive-frequency  field  operator  a{u!,z),  represent¬ 
ing  entangled  photons,  evolves  inside  the  crystal  and  may  be  written  in  terms  of  an  envelope 
b{u),z)  which  varies  slowly  along  the  propagation  direction:  a{u,z)  —  b{uj,  z)  exp[ik{ui)z]. 
The  propagation  equations  describing  entangled  photon  generation  are 
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db{u^  z)  _  [2* A/:(a;)^] , 

dz 


db^u,  z) 
dz 


=  —iK*{uj)b{ujp  —  cj,  z)  exp  [—iAk{(jj)z]. 


(7.1) 


where  k{u})  and  Ak{uj)  are  the  coupling  factor  and  wave-vector  mismatch,  respectively.  The 
solution  for  the  output  field  at  ^  =  Z/,  expressed  in  terms  of  the  vacuum  field  avac(^)  at  the 
input  of  the  crystal,  is 


aout(^^)  =  A{u))ay^{Lo)  +  -  a;),  (7.2) 

where,  to  preserve  the  commutation  relations,  the  functions  A{u;)  and  B{u;)  satisfy  |A(a;)p  — 
\B{u;)\'^  =  1  and  A{iv)B{iOp  —  to)  =  B{to)A{ijp  —  u). 

The  time-domain  output  field  oi)erator  is  related  to  its  frequency-domain  counterpart 
[Eq.  (7.2)]  by  the  inverse  Fourier  transform,  aout(0  =  «out(^)  exp(— za;t)dc<;,  and  is  nor¬ 

malized  so  that  the  total  rate  of  generated  photons  exiting  the  crystal  is  i?out  =  (^lut(^)^out(^))- 
The  generated  photons  are  separated  into  two  channels,  denoted  as  Channel  1  and  Chan¬ 
nel  2,  using  a  50/50  beam  splitter.  The  field  operators  at  the  outputs  of  the  beam  splitter 
are  ai(t)  =  a2{t)  =  -^aoutit)-  The  photons  are  modulated  by  periodic  phase  modula¬ 
tors  whose  time-domain,  Fourier-series  transfer  functions  are  mi{t)  =  Qk  ^^p{—ikuJmi)  in 
Channel  1  and  m2{t)  =  r/ exp(— in  Channel  2,  with  Fourier  transforms  mi  (a;)  = 
—  kujm)  and  m2(cu)  =  respectively.  With  the  *  symbol  denot¬ 

ing  convolution,  the  frequency-domain  modulated  fields  are  ai(a;)  =  ai{uj)  *  mi  (a;)  and 
n2(a^)  =  fi2{uj)  *  7712(00).  Substituting  ai(u;),  a2{oo),  77ii{oo),  and  7712(00)  into  the  expressions 
for  Si  (a;)  and  0.2(00)  yields 


Si  (a;) 


d2(co) 


1 

^  ^  [-^(^  ktOfn)Oy;^c(cO  kcOm) 

^  k=—oo 

+  B(lO  kLOjji)o\^(LOp  +  /lCi^jti)]  , 

— ^  ^  ^  Tl  [.4(cJ  ^^m)nvac(a^  ^^m) 

^  /=— 00 

“h  B(uJ  lLOfYi)o^Qj^(uJp  +  ^^m)]  • 


(7.3) 
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The  modulated  photons  are  frequency  correlated  by  passing  each  through  identical  monochro¬ 
mators  whose  output  slits  may  be  translated  to  select  frequencies  uji  =  (3x\  in  Chan¬ 
nel  1  and  UJ2  =  0^2  in  Channel  2,  where  the  constant  /?  is  the  linear  dispersion  of  the 
grating  systems.  The  monochromators  (spectral  filters)  have  field  transmission  functions 
Hi{uj  —  Pxi)  and  H2{u)  —  (3x2)-  The  filtered  field  operators  in  Channels  1  and  2  are 
dif(o;, xi)  =  —  Pxi)  and  a2f{uj,X2)  =  a2{uj)H2{u>  —  0X2),  respectively.  The  count 

rates  at  the  outputs  of  the  monochromators  are  given  by  Ri{xi)  =  2:i)aif(t, 2:1))  and 

R2{x2)  =  {al({t,X2)a2{{t,X2)).  These  rates  are 


1  poo 

Ri{xi)  =  —  '^\Qk\'^  \B{u;-ku;^)f\Hi{u;-Pxi)fdu, 


1  P^ 

R2{x2)  =  J-  /  \B{U- 


\H2{oJ  -  0X2)fdoj. 


(7.4) 


Assuming  a  gate  width  T,  the  coincidence  rate  for  the  two  detectors  is  related  to  the 
second-order  Glauber  correlation  function  Xi,t2,X2)  = 

{a\^{t2,X2)o\^{ti^xi)aif{ti^Xi)a2f{t2^X2)).  With  the  assumption  that  the  resolution  of  the 
monochromators  is  high,  or  equivalently  that  the  filter  widths  are  small  (as  compared  to  the 
modulation  frequency  o;^),  it  can  be  shown  that  the  correlation  function  depends  only  on 
the  difference  of  the  arrival  times  r  =  t2  —  t\,  and  the  coincidence  rate  is 


rT/2 

Rc{xi,X2)=  /  G^‘^\r,Xi,X2)dr. 

J-T/2 


-T/2 

Equation  (7.5)  may  be  expanded  using  Wick’s  theorem  and  shown  to  be  given  by 


(7.5) 


Rc{XuX2)  =  Ri{Xi)R2{x2)T  +  / 

J  —c 


'^qkrn-kFkix,  Xi,X2) 


2 

dr, 


k=—(X) 


(7.6) 


where  A  =  0{xi  +  X2)  —  cOp,  n  =  [A/uz^  +  and 


1  /■°° 

Fk{T,  Xi,X2)  =  1-  /  A{u;  -  ku;rn)B{u;p-u;  +  kum) 

47r  Jo 

X  Hi{u!  —  0xi)H2{u!p  —  uj—  0X2  +  TWm)  exp{iu!T)du}.  (7.7) 


35 


The  first  term  in  Eq.  (7.6)  is  the  result  of  accidental  coincidences  between  unpaired  pho¬ 
tons  in  a  gate  width  T.  The  second  term  is  the  coincidence  rate  between  paired  photons 
and  captures  the  modulation  effects  described  in  this  section.  To  obtain  Eqs.  (7.4)-(7.7),  we 
have  assumed  that  the  transmission  widths  of  the  monochromators  are  small  as  compared  to 
the  modulation  frequency  and  large  as  compared  to  the  inverse  of  the  temporal  gate  width 
T.  In  our  experiment  these  assumptions  are  satisfied  by  factors  of  3.5  and  11,  respectively. 

If  we  further  assume  that  A{uj)  and  B{ujp  —  uj)  are  constant  in  the  vicinity  (±150  GHz  in 
Fig.  15)  of  a;  =  I3x\  and  are  equal  to  Aq  and  jBq,  respectively,  then  Eq.  (7.6)  becomes 


—  R1R2T  ±  CnH^‘^\nuJrn  ~  ^)) 
where  =  \Hi{uj)\‘^  *  |i^2(t^)P,  and 


Ri 

R2 

Cfi 
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47r 
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Stt 


/oo 

■00 

/OO 

-OO 

OO  ^ 

^0^0  ^  ^  Qk'^n—k  • 


(7.8a) 

(7.8b) 

(7.8c) 


The  solid  curves  in  Fig.  15  are  theoretical  fits  to  the  data  using  Eq.  (7.8)  shifted  horizon¬ 
tally  so  as  to  match  center.  The  Fourier  series  coefficients  for  sinusoidal  phase  modulators 
are  Bessel  functions  with  qk  =  Jk{—^i)  and  r/  =  J/(— ^^2),  where  Si  and  62  are  the  modu¬ 
lation  depths  in  Channels  1  and  2,  respectively  {\Si\  =  |(52|  =  1.5  in  our  experiment).  We 
model  the  monochromator  response  functions  in  Channels  1  and  2  as  Gaussians  with  FWHM 
bandwidths  F:  Hi{uj)  =  ai  exp  [— 2  ln(2)a;^/r^]  and  H2{uj)  =  a2  exp  [— 2  ln(2)c<;^/r^].  (The 
monochromator  in  Channel  1  is  the  mirror  image  of  the  one  in  Channel  2  which  has  a  mea¬ 
sured  FWHM  bandwidth  of  8.5  GHz.)  The  transfer  functions  include  fitting  parameters  ai 
and  a2  used  in  Fig.  15  to  account  for  transmission  losses  and  the  difference  in  detection 
efficiencies  of  the  photon  counters. 


36 


To  obtain  the  constants  Aq  and  Bq,  for  each  case  in  Fig.  15,  we  measure  the  average 
value  of  and  use  Eq.  (7.8b)  to  calculate  |i?o|-  We  obtain  |j4o|  from  the  commutator¬ 
preserving  condition  \Aq\^  —  |i?oP  =  !•  For  all  curves,  the  fitting  parameters  are  taken  as 
a\  =  1.20  X  10~^  and  =  5.59  x  10~^.  These  values  are  in  good  agreement  with  loss 
measurements  and  estimates  of  the  photon  counter  detection  efficiency,  where  we  note  that 
the  id400  detector  in  Channel  1  has  a  detection  efficiency  an  order  of  magnitude  larger  than 
the  SPCM-AQR-16-FC  detector  in  Channel  2. 


8  Chirp  and  Compress  with  Biphotons 

Background 

When  the  bandwidth  of  an  optical  pulse  is  equal  to,  or  greater  than,  its  central  frequency 
and  when  all  of  its  spectral  components  are  in  phase,  the  pulse  develops  a  characteristic 
waveform  that  is  termed  as  single  cycle.  A  chirped,  quasi  phase  matched  nonlinear  crystal 
is  a  technique  for  generating  nonclassical  pairs  of  photons  (biphotons)  whose  characteristic 
coincidence  time,  as  measured  at  distant  detectors,  is  a  single  optical  cycle.  The  elements  of 
the  technique  are  (1)  the  suggestion  for  using  parametric  down  conversion  to  spontaneously 
generate  pairs  of  entangled  photons  whose  instantaneous  frequencies  are  chirped  in  opposite 
directions,  and  (2)  the  use  of  the  nonlocal  nature  of  entangled  photons  to  allow  the  dispersion, 
as  experienced  by  one  photon,  to  cancel  out  the  dispersion  of  the  second  photon  and  to 
compress  the  biphoton  wave  packet. 

A  key  motivation  for  the  study  of  single-cycle  biphotons  is  their  potential  application  to 
nonlinear  optical  processes  with  nonclassical  fields.  One  example  is  the  use  of  sum  frequency 
generation  (SFG)  as  an  ultrafast  correlator  of  biphotons.  The  efficiency  for  generating  sum 
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frequency  photons  varies  inversely  with  the  width  of  the  incoming  biphoton;  i.e.,  single-cycle 
biphotons  behave  as  if  they  have  an  effective  power  equal  to  their  energy  divided  by  their 
temporal  width.  Similarly,  in  the  absence  of  intermediate  resonances,  single-cycle  biphotons 
maximize  the  two  photon  transition  probability.  Other  uses  for  ultra-wideband  biphotons 
include  nonclassical  metrology  and  large  bandwidth  quantum  information  processing. 


Figure  16:  Chirped-poled  nonlinear  crystal.  Arrows  indicate  the  sign  of  the  nonlinear  coef¬ 
ficient. 

As  shown  in  Fig.  16,  we  make  use  of  a  quasi-phasematched  (QPM)  periodically  poled 
crystal.  The  up-down  arrows  show  the  sign  of  the  nonlinear  coefficient  in  each  domain  of 
the  crystal.  These  domains  are  reversed  with  a  period  such  that  the  corresponding  spatial 
frequency,  27r/A,  is  linearly  chirped.  The  poling  period  is  chosen  so  that  the  signal  frequency 
is  phase  matched  for  red  emission  at  the  left  end  of  the  crystal  and  for  blue  emission  at  the 
right  end  of  the  crystal.  Paired  photons  that  are  emitted  from  the  right  end  arrive  at  their 
respective  photodetectors  at  the  same  time.  Paired  photons  that  are  emitted  from  the  left 
end  of  the  crystal  arrive  at  the  photodetector  with  a  time  difference  determined  by  their 
group  velocities.  In  the  ideal  case  where  there  is  no  group  velocity  dispersion,  the  biphoton 
wave  packet  is  (exactly)  linearly  chirped.  Figure  17  shows  the  near-white  light  spectrum  of 
a  2  cm  long  QPM  crystal  of  LiNbOa  pumped  with  a  monochromatic  laser  at  a  wavelength  of 
0.420  iim.  Here,  the  spatial  frequency  of  the  domain  reversals  varies  linearly  with  distance 
and  is  chosen  so  that  the  crystal  is  phase  matched  at  a  signal  wavelength  of  0.750  //m  at  the 
left  end  and  0.464  /xm  at  the  right  end.  This  corresponds  to  a  poling  period  of  3.11  fim  at 
the  left  end  and  7.02  //m  at  the  right  end  of  the  crystal.  When  compressed,  this  spectrum 
corresponds  to  a  biphoton  with  a  temporal  length  that  is  nearly  a  single  optical  cycle  at  the 
degenerate  wavelength  of  0.840  /xm.  The  refractive  index  as  a  function  of  frequency,  for  this 
figure,  is  ne{(jj)  and  is  obtained  from  the  Sellmeier  equation  for  e-polarized  light  in  LiNbOa. 
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Figure  17:  Calculated  spectrum  of  spontaneously  generated  chirped  photons  for  a  2  cm  long 
QPM  crystal  of  LiNbOa  pumped  at  0.420  //m. 

Theoretical  Work 

We  have  developed  a  strong  understanding  of  the  physics  of  nonlinear  optics  with  biphotons, 
specifically  sum  frequency  generation,  and  now  have  a  well-understood  theory  to  solve  other 
linear  and  nonlinear  quantum  optical  processes  using  down-converted  photons.  We  have 
developed  equations  and  numerical  codes  which  allow  us  to  simulate  down  conversion  in 
periodically-  and  chirped-poled  nonlinear  crystals.  Of  particular  importance  are  codes  which 
simulate  biphoton  generation  in  poled  crystals  with  arbitrary  domain  patterns,  leading  to 
new  ideas  for  obtaining  very  short  biphotons  (e.g.  using  a  nonlinearly-chirped  crystal  to 
compensate  for  dispersion).  Our  numerical  tools  also  demonstrate  the  effects  of  the  rather 
large  (and  unavoidable)  poling  errors  present  in  real  crystals  and  have  lead  to  a  particularly 
interesting  discovery  that  the  chirp/compress  strategy  for  generating  single-cycle  biphotons 
is  very  robust  to  these  errors.  As  an  example  showing  the  usefulness  of  our  code,  we  apply 
it  to  an  experiment  in  progress.  We  are  currently  setting  up  an  ultrafast  correlator  to 
measure  the  width  of  a  modestly-chirped  biphotoii.  We  are  using  a  chirped-poled  crystal 
from  HC  Photonics  which  has  a  linear  k- vector  chirp.  The  crystal  is  magnesium  oxide-doped 
stoichiometric  lithium  tantalate  (MgO:SLT),  with  a  poling  period  which  varies  from  8.0481 
to  8.0223  /jtm  over  its  2  cm  length  from  input  to  output.  Our  numerical  code  computes  the 
theoretical  biphoton  wavefunction  generated  by  the  crystal,  as  shown  in  Fig.  18. 
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Figure  18:  Numerical  simulation  computing  the  output  of  a  chirped  nonlinear  crystal  (2  cm 
long  MgOiSLT).  Calculated  spectrum  of  the  signal  (left)  and  idler  (right)  photons. 

As  a  check  to  our  theory,  we  have  measured  the  spectra  of  the  signal  and  idler  photons 
generated  by  the  MgOiSLT  crystal,  shown  in  Fig.  19.  The  measurements  were  made  with  a 
Jarrell  Ash  MonoSpec  27  CCD  spectrograph.  The  bandwidth  appears  to  be  accurate,  and 
some  of  the  features  (e.g.  lobes)  predicted  by  the  numerical  code  can  be  seen  in  the  spectra. 

Measured  signal  spectaim  Measured  idler  spectrum 


Figure  19:  Measured  spectra  of  the  signal  and  idler  photons  generated  by  a  chirped  nonlinear 
crystal  (2  cm  long  MgO'.SLT). 

To  measure  the  correlation  function  for  the  biphotons  generated  in  the  chirped  crystal, 
we  plan  to  use  a  1  mm  LiNbOs  crystal  to  generate  532  nm  SFG  from  the  biphotons  as  we 
vary  the  delay  between  them  (more  details  will  be  discussed  in  the  experimental  section). 
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The  numerical  code  also  predicts  the  SFG  correlation  curves  for  such  an  experiment,  which 
are  shown  in  Fig.  20  for  the  chirped  (left)  and  compressed  (right)  biphotons.  The  simulation 
which  produced  Figs.  18  and  20  did  not  include  poling  errors,  but  they  are  easily  simulated 
with  the  same  code. 

Quantum  SFG  rate  Quantum  SFG  rate 

(uncompressed,  =  1 .003  (mm))  (compressed,  =  1 .003  (mm)) 


Figure  20:  Numerical  calculation  of  expected  SFG  correlation  measurement.  Left:  corre¬ 
lation  function  of  the  chirped  biphoton.  Right:  correlation  function  for  the  biphoton  after 
compression. 

In  addition  to  the  theoretical  work  directly  related  to  our  chirp  and  compress  experiment, 
we  have  also  studied  other  problems  related  to  broadband  biphotons.  We  have  introduced  a 
new  form  of  spectroscopy,  called  Quantum  Fourier  Transform  Spectroscopy  (QFTS),  which 
allows  determination  of  the  complex  refractive  index  of  an  unknown  material  using  a  nar¬ 
rowband  detector.  The  technique  relies  on  sum  frequency  generation  of  biphotons  and  ho- 
modyning  the  coherent  output  with  a  local  oscillator;  it  cannot  be  done  with  classical  light 
sources.  By  working  in  the  Heisenberg  picture,  our  theory  has  made  it  especially  straight¬ 
forward  to  simulate  parametric  down  conversion  in  the  high  gain  regime,  allowing  the  same 
simple  equations  to  be  used  to  simulate  nonlinear  and  quantum  optics  with  the  N-photon 
states  generated  in  this  regime.  One  interesting  result  of  this  analysis  is  that  nonlocal  can¬ 
cellation  of  dispersion  (a  property  which  allows  biphotons  to  be  compressed  by  acting  on 
only  one  photon  of  the  pair)  is  observable  at  classically  high  power  levels.  Our  theory  has 
also  led  us  to  predict  Hong-Ou-Mandel  interference  with  chirped  biphotons.  In  the  case  of 
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ideal  linear  chirps  (no  matter  how  large),  the  shortest  possible  HOM  dip  can  be  achieved 
with  zero  dispersion  compensation.  This  makes  broadband  biphotons  generated  by  chirped 
crystals  attractive  for  experiments  which  utilize  narrow  HOM  dips  (such  as  quantum  optical 
coherence  tomography). 

Experimental  Work 


Idler 


Figure  21:  Sum  frequency  generation  correlator  setup  using  a  dichroic  mirror  to  separate 
signal  and  idler  photons.  The  biphoton  is  compressed  in  the  signal  channel,  and  the  idler 
photon  is  delayed  for  correlation. 

We  are  currently  in  the  process  of  setting  up  a  proof-of-principle  experiment  to  demon¬ 
strate  the  idea  of  chirping  and  subsequently  compressing  biphotons  to  achieve  narrow  cor¬ 
relation  widths  and  increased  efficiency  of  a  nonlinear  process.  The  basic  correlation  setup 
is  shown  in  Fig.  21.  The  biphotons,  incident  from  the  left,  separate  into  signal  and  idler 
photons  at  a  dichroic  beam  splitter.  The  beam  splitter  is  a  long-pass  edge  filter  (Semrock 
LP02-1064RS-25)  with  a  very  sharp  transition  edge  that  can  be  angle  tuned  to  1064  nm,  the 
wavelength  separating  our  photons  (see  Fig.  18).  The  idler  photon  passes  through  the  filter 
and  is  retroreflected  back  through  the  filter.  The  retroreflector  is  mounted  on  an  automated 
translation  stage  which  can  travel  up  to  24  mm  in  10  nm  steps  (Newport  VP-25XL),  allowing 
high-resolution  correlation  traces.  The  signal  photon  is  reflected  from  the  filter  and  passes 
through  a  separate  arm  in  the  correlator.  In  the  first  phase  of  the  experiment,  this  arm 
will  simply  return  the  photon  to  be  overlapped  with  the  delayed  idler  photon  at  the  beam 
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splitter.  In  the  second  phase  of  the  experiment,  a  compressor  (custom  chirped  mirrors  from 
Femtolasers)  will  be  inserted  in  the  signal  beam  path  to  compress  the  biphoton  chirp.  The 
recombined  signal  and  (delayed)  idler  photon  are  focused  into  a  1  mm  long  periodically-poled 
MgO-doped  LiNbOa  crystal  to  generate  532  nm  SFG.  The  SFG  will  be  measured  by  a  single 
photon  counting  module  (PerkinElmer  SPCM-AQR-16-FC)  and  recorded  as  a  function  of 
idler  photon  delay  to  measure  the  correlation  function  of  the  biphoton.  When  the  two  phases 
of  the  experiment  are  complete,  we  expect  to  see  compression  of  the  correlation  width  by  a 
factor  of  10,  with  a  corresponding  increase  in  SFG  rate  (as  predicted  by  Fig.  20). 


9  Nonlinear  Optics  at  X-Ray  Wavelengths 

Our  motivation  for  this  work  is  the  construction  of  the  new  LCLS  (Linac  Coherent  Light 
Source)  at  Stanford.  This  source  will  use  about  1/3  of  the  present  SLAG  accelerator  and 
will  produce  PEL  radiation  at  wavelengths  as  short  as  1.5  Angstroms.  The  projected  op¬ 
erating  parameters  for  this  light  source  are:  average  power=  lO^'^  photons  per  second; 
pulse  duration^  2  fs;  photons  per  pulse=l.l  x  10^^;  repetition  rate  =100  Hz;  angular 
divergence=0.8  x  10”^  radians.  Taken  together  these  parameters  should  allow  the  first 
demonstration  of  new  types  of  nonlinear  x-ray  effects. 

We  anticipate  LCLS  experiments  in  the  following  areas:  1)  The  demonstration  of  second 
harmonic  generation  in  its  own  right  and  for  application  to  the  measurement  of  the  temporal 
structure  of  the  LCLS  femtosecond  time  scale  laser  pulses.  As  in  the  visible  region  of  the 
spectrum  these  experiments  will  be  done  by  splitting  the  x-ray  beam  and  temporally  delaying 
one  portion  relative  to  the  other,  thereby  allowing  the  use  of  a  slow  detector.  2)  The 
demonstration  of  mixing  of  a  fs  time  scale  visible  laser  with  the  LCLS  beam.  We  have  found 
that  the  conversion  efficiency  for  this  latter  process  is  about  four  orders  of  magnitude  larger 
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than  for  x-ray  frequency  doubling,  and  it  is  therefore  likely  that  this  effect  will  allow  the 
development  of  practical  correlation  with  a  resolution  of  the  shortest  available  Ti:Sapphire 
lasers.  3)  Because  this  type  of  non-linear  x-ray- visible  process  is  dependent  on  the  locations  of 
the  valence  electrons  it  is  possible  that  it  can  be  developed  into  a  method  for  the  studying  the 
shape  of  chemical  bonds.  4)  Experiments  on  parametric  down  conversion:  These  experiments 
would  first  be  aimed  at  studying  the  basic  process,  for  example  the  correlated  count  rate  as 
a  function  of  angle,  position,  and  frequency.  5)  More  advanced  experiments  will  be  aimed 
at  measuring  the  temporal  length  of  the  generated  correlated  biphotons.  Of  importance 
our  calculations  show  that  these  photons  are  born  with  a  length  of  less  than  10“^^  sec.  It 
may  therefore  be  surprising  easy  to  generate  photons  allowing  measurement  of  atomic  scale 
distances.  6)  It  is  likely  that  the  LCLS  will  lead  to  a  variety  of  experiments  which  make  use 
of  entanglement  at  x-ray  wavelengths.  These  include  studies  of  effects  such  as  ghost  imaging 
and  Hong-Ou-Mandel  beam  splitter  experiments  on  ultrashort  x-ray  correlation 

In  the  following  we  summarize  our  theoretical  study  and  describe  our  proposed  experi¬ 
ments. 


X-Ray  Nonlinear ites 

The  central  concept  pointed  out  in  the  earlier  papers  of  Freund  and  Levine  {Phys.  Rev,  Lett. 
23,  854,  1969)  and  Eisenberger  {Phys,  Rev.  Lett.  26,  684,  1971)is  that,  since  X-ray  photons 
have  energies  that  are  large  as  compared  to  the  electron  binding  energy  of  light  elements,  but 
still  small  as  compared  to  the  electron  rest  mass,  that  the  nonlinearity  of  an  element  such 
as  diamond  may  be  calculated  by  treating  all  of  the  electrons  in  the  atom  as  free  particles; 
and  therefore  representing  the  nonlinear  medium  as  a  very  dense,  (on  order  of  10^^  electrons 
per  cubic  centimeter)  cold  plasma.  Of  most  importance,  the  periodic  lattice  structure  causes 
the  plasma  to  be  periodic.  In  a  homogeneous  plasma  the  projection  of  the  nonlinear  current 
density  against  the  driving  fields  is  zero  and  the  effective  second  order  nonlinearity  is  zero. 
The  charge  density  of  this  effective  plasma  is  determined  by  the  reciprocal  lattice  k- vectors, 
weighted  by  the  appropriate  structure  and  form  factors.  (Note  that  this  is  different  from 
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conventional  nonlinear  optics  in  a  periodically  poled  medium  where  the  nonlinearity  exists 
irrespective  of  the  poling  period.) 

The  nonlinearity  of  this  dense  plasma  is  determined  by,  at  first,  seemingly  different  phys¬ 
ical  processes.  These  are  a)  the  Lorentz-force,  b)  convection,  as  determined  by  the  spatial 
variation  of  velocity,  and  c)  a  current  density  caused  by  the  spatial  variation  of  the  time- 
varying  charge  times  the  first  order  velocity.  In  recent  work,  we  have  shown  that  the  plasma 
nonlinearity  will  satisfy  detailed  balance  and  the  Manley  Rowe  relations.  This  is  at  first 
not-obvious;  for  example  one  may  show  that  the  Lorentz-force  term,  when  taken  separately 
from  the  other  terms  does  not  satisfy  the  Manley- Rowe  relation  and  would  not,  in  its  own 
right  produce  correlated  photons.  But  by  taking  the  sum  of  the  Lorentz  term,  the  convec¬ 
tion  term,  and  the  second  order  current  term,  all  on  an  equal  footing,  detailed  balance  is 
regained. 

It  should  perhaps  be  noted  that  though  the  nonlinear  coefficients  at  x-ray  frequencies  are 
much  smaller  than  those  at  optical  wavelengths,  this  is  offset  by  the  the  ability  to  focus  the 
driving  beams  much  more  tightly.  Also,  in  the  case  of  parametric  down  conversion,  though 
the  conversion  efficiency  from  the  vacuum  field  is  reduced,  the  number  of  modes  in  k-space 
is  greatly  increased. 

We  noted  that  the  second  order  nonlinear  effects  we  discuss  here  can  be  considered  as  as 
a  combination  of  an  electric  quadrupole-type  and  a  magnetic  dipole-type  nonlinearities.  As 
was  shown  by  P.  S.  Pershan  (Phys.  Rev.  130,  919  (1963))  these  type  of  nonlinearities  are  not 
zero  in  materials  having  a  center  of  inversion;  i.e.  second  order  nonlinear  x-  ray  processes 
may  occur  in  materials  possessing  inversion  symmetry. 


SHG  at  X-Ray  frequencies 

9.0.1  N  onlinear  ity 

For  SHG  process  the  nonlinear  current  density  is  (Y.  R.  Shen  The  Principles  of  Nonlinear 
Optics  p.  543  ) 
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Where  q,m  are  the  electron  charge  and  mass  respectively,  uj  is  the  frequency  of  the  funda¬ 
mental  beam,  po  is  electron  chrage  density  (in  the  x-ray  regime  this  is  the  Fourier  component 
corresponding  to  the  lattice  reciprocal  vector),  is  the  electric  field  of  the  the  fundamental 
beam,  and  is  the  direction  of  the  electric  field  of  the  second  harmonic  beam. 

9.0.2  Phase  Matching  for  a  Single  Pump  Beam 

The  phase  matching  diagram  for  the  Lane  geometry  is  shown  in  Fig.  22  With  and  fc2a;5and 
Q  denoting  the  k-vectors  of  the  fundamental,  second  harmonic,  and  reciprocal  lattice  vector, 
the  phase  matching  condition  for  a  second  harmonic  generation  (SHG)  process  at  x-ray 
wavelengths  is  2k^  +  Q  =  k2uj- 


Figure  22:  Phase  matching  diagram  for  x-ray  SHG  for  a  single  fundamental  beam,  k^,  k2uj, 
are  the  wave  vectors  of  the  fundamental  and  second  harmonic  beams  respectively.  Q  is  the 
reciprocal  lattice  vector  of  the  diamond  crystal 

As  is  shown  in  Fig.  22,  we  define  and  O^uj  as  the  angles  between  the  atomic  planes 
and  the  fundamental  beam  and  the  second  harmonic  beam  respectively.  Ignoring  dispersion, 
2k^  =  k2u}  and  the  phase  matching  condition  is  identical  to  that  of  Bragg’s  law  for  an  x-ray 
beam  at  a  frequency  of  2(jj.  However,  dispersion  does  exist  and  02a;  and  02a;  have  values 
which  are  close  to  but  different  from  the  Bragg  angle  (at  02a;).  Of  importance  the  SHG  and 
Laue  rocking  curves  do  not  overlap  ,  thereby  avoiding  diffraction  of  the  generated  second 
harmonic  beam.  For  the  (111)  reflection  of  diamond  and  for  a  second  harmonic  signal  at  16 
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keV,  the  deviation  of  6201  from  the  Bragg  angle  is  24  /x  rad.  Using  the  dynamical  theory  of 
x-ray  diffraction  in  the  plane  wave  approximation  (Dynamical  Theory  of  X-Ray  Diffraction, 
Andre’  Authier,  ’’Oxford  Science  Publication”,  New  York,  2001),  one  finds  that  the  Darwin 
width  at  16  keV  is  12  /x  rad.  Recently,  by  solving  numerically  the  full  coupled  equations 
(including  SHG  and  linear  Lane  diffraction),  we  found  that  for  a  fundamental  spot  size  of 
2//m  and  assuming  phase  matching  for  the  SHG  ,  less  than  one  percent  of  the  generated 
second  harmonic  is  linearly  diffracted.  We  therefore  neglect  the  linear  diffraction  of  the 
second  harmonic  beam. 


9.0.3  Beam  Profiles  and  Efficiency 

By  using  the  undepleted  pump  approximation  and  neglecting  the  scattering  of  the  generated 
second  harmonic  signal,  together  with  the  nonlinear  susceptibility  calculations,  we  have  car¬ 
ried  out  calculations  that  predict  the  SHG  efficiency  and  its  dependance  on  the  propagation 
distance.  We  also  calculated  the  intensity  profile  in  the  transverse  plane  of  the  generated 
second  harmonic  beam.  The  slowly  varying  envelope  equation  at  the  second  harmonic  in  the 
Laue  configuration  and  for  coplanar  electric  fields  reads  as 


27r 


sin  9  b 


dx 


+  cos  9  b 


dz  ) 


i  (x'(2w)  -  x'{^)  cos{29b)  +  2  cos{29 b)S9i  +  ix"(2cj)^  E^^^ 

(9.2) 


poq^  sin(20B)  [  4cos(20b)-1  1 

Here  ^  is  the  nonlinear  susceptibility, 0b,  is  the  vacuum 

Brag  angle  at  the  second  harmonic,  A  is  the  vacuum  wavelength  at  the  second  harmonic, 
X'(2u;),x\^)  ^re  the  real  parts  of  linear  susceptibilities  of  the  second  harmonic  and  the  fun¬ 
damental  respectively.x^^(2ct;)  is  the  imaginary  part  (loss)  of  the  linear  susceptibility  of  the 
second  harmonic  beam.  S0i  is  the  angular  deviation  of  the  pump  beam  from  the  vacuum 
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Figure  23:  (a) Normalized  Intensity  of  the  fundamental  vs.  x  and  z  (b)  SHG  efficiency  vs.  x 
and  z 

Bragg  angle.  The  absorption  of  the  pump  is  introduced  by  takin  the  pump  field  as  a  decaying 
exponent. 

The  efficiency  of  the  SHG  process  depends  linearly  on  the  intensity  of  the  fundamental 
beam.  Therefore  we  will  focus  the  LCLS  beam  to  a  beam  spot  of  a  few  microns  .  Though  this 
yields  the  highest  efficiency,  it  causes  the  second  harmonic  beam  to  be  generated  at  an  angle 
of  Ou;  +  O2U}  with  regard  to  the  fundamental  driving  beam.  As  a  result,  the  beams  overlap  in 
space  for  a  distance  which  is  determined  by  the  spot  size  of  the  fundamental  beam  and  the 
angle  between  the  two  beams.  For  a  pumping  beam  diameter  of  2/im  at  8  Kev,  this  distance, 
which  we  will  refer  as  the  “interaction  length”,  is  10  /im.  For  transparent  crystals  which 
are  thicker  than  interaction  length,  the  SHG  efficiency  grows  as  the  propagation  length,  and 
not  as  its  square,  as  in  conventional  nonlinear  optics  (or  in  the  case  where  this  walk-off  effect 
is  negligible).  This  effect  results  in  a  reduction  of  the  power  of  the  the  generated  second 
harmonic  and  a  spreading  of  the  second  harmonic  signal  in  the  transverse  plane.  This  effect 
is  shown  in  Fig.  23.  Fig.  23(a)  shows  the  normalized  intensity  of  the  fundamental  vs.  x  (a 
coordinate  parallel  to  the  crystal’s  boundary)  and  z  (  a  coordinate  normal  to  the  boundary). 
Fig.  23(b)  shown  the  efficiency  of  the  SHG  vs.  x,  z  in  normalized  units. 

To  estimate  the  efficiency  we  assume  a  fundamental  flux  of  4  x  10^  photons  per  pulse  at 
8keV,  a  pulse  length  of  2fs  and  beam  diameter  of  2  microns.  We  have  used  the  diamond 
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z(|Lim) 


Figure  24:  SHG  efficiency  vs.  crystal  thickness. 

(Ill)  reflection  for  phase  matching  the  nonlinear  medium.  Fig.  25  shows  the  efficiency  as 
function  of  the  diamond  thickness.  For  short  crystal  thicknesses  (  but  thicker  than  the  inter¬ 
action  length)  the  efficiency  increases  linearly  with  the  thickness,  but  at  larger  propagation 
distances,  absorption  plays  a  limiting  important  role,  and  the  efficiency  dependance  on  the 
propagation  distance  becomes  sub-linear.  The  intensity  profile  for  several  crystal  thickness 
are  shown  in  Fig  24.  One  can  see  very  clearly  that  the  intensity  profile  is  broaden  with  the 
crystal  thickness.  While  for  thickness  of  20fim  the  beam  diameter  is  still  comparable  to  the 
beam  diameter  of  the  fundamental,  for  a  thickness  of  lOO/im  the  profile  is  5  times  broader, 
and  for  a  500/xm  the  profile  is  about  100  times  broader  than  the  intensity  profile  of  the  pump 
beam. 


9.0.4  Experimental  Setup 

The  anticipated  second  harmonic  data  rate  is  on  the  order  of  a  single  to  ten  photons  per  pulse 
when  LCLS  is  operating  in  short  pulse  mode  and  the  beam  is  focused  to  a  2  micron  spot  size. 
The  expected  signal  is  extremely  small  and  background  levels  must  be  carefully  considered. 
As  such,  a  crystal  analyzer  will  be  used  to  isolate  the  generated  second  harmonic  photons 


49 


Figure  25:  Normalized  profile  of  the  SHG  efficiency  for  several  crystal  thickness  (a)  20/xm 
(b)  100/im,  (c)  250/txm,  and  (d)  500/im 

from  the  elastic  and  inelastic  background  signal  from  fundamental  and  3^^  harmonic  photons 
in  the  incident  beam.  A  channel  cut  analyzer  is  needed  to  prevent  background  signal  from 
double  scattering  events  from  the  sample  and  first  analyzer  crystal.  However,  this  scheme  is 
dependent  upon  an  incident  beam  that  is  free  from  2^^  harmonic  photons.  This  is  achieved 
using  a  combination  of  a  silicon  111  monochromator,  harmonic  rejection  mirror  system  and 
the  chromaticity  of  the  XPP  focusing  lens  system.  This  combination  of  X-ray  optics  will 
create  a  spectrally  pure  beam  of  sufficient  quality  for  the  proposed  experiment.  A  schematic 
of  the  setup,  including  photon  numbers  after  each  optical  element,  is  displayed  in  upper  part 
of  Figure  26. 
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2  fs,  20  pC  LCLS  Operation  at  8  keV 


Figure  26:  Upper:  Schematic  of  the  SHG  experimental  setup.  A  silicon  111  channel  cut 
monochromator,  2  bounce  harmonic  rejection  mirror  system  and  beryllium  compound  re- 
fractive  focusing  lens  will  be  used  to  prepare  the  incident  X-ray  beam.  The  expected  photon 
numbers  for  the  fundamental,  2^^  and  3^^  harmonic  are  displayed  before  the  monochromator, 
after  the  monochromator  and  after  the  Harmonic  Rejection  Mirror  system.  A  silicon  111 
analyzer  crystal  and  photodiode  will  be  used,  with  apertures,  to  detect  the  generated  second 
harmonic  signal.  Lower:  Schematic  diagram  for  x-ray  intensity  correlation. 
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9.0.5  Correlation  of  the  X-Ray  Free  Electron  Laser  (XFEL) 


The  LCLS  is  expected  to  deliver  femtosecond  pulses  with  the  high  powers  that  are  sufficiently 
high  to  allow  observation  of  SHG.  In  the  optical  regime,  SHG  is  now  widely  used  in  auto¬ 
correlators  for  measuring  the  pulse  duration  of  sub  100  femtosecond  pulses.  It  is  therefore 
reasonable  to  believe  that  a  similar  technique  may  be  established  in  the  x-ray  region.  Lower 
part  of  Fig.  26  shows  a  possible  scheme  for  autocorrelation  measurement.  The  beam  from 
the  LCLS  will  be  split  into  two  portions. One  portion  of  the  beam  will  propagate  through 
a  delay  line  and  then  the  two  portions  of  beam  will  recombined  at  the  nonlinear  material. 
A  second  harmonic  beam,  proportional  to  the  product  of  the  intensities  of  two  portion  of 
the  fundamental  beam  will  be  generated.  The  slow  detector  will  integrate  the  intensity  with 
respect  to  the  time  .  By  varying  the  delay  time  between  the  two  portion  of  the  beam,  one 
can  obtain  the  autocorrelation  function  of  the  beam  intensity. 
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B  Chirped  Crystal  Simulation  Code 


clear  ell; 
tic; 

foraat  long; 

eet(0, 'def aultazeef onteize' , 12) ; 
eet(0, 'defeultteztf onteize' , 12) ; 

%  Conet ante 
c  -  3e6; 

hb  -  6.83e-34/(2*pl); 
epeO  -  le-9/(36*pl) ; 
eteO  -  377; 

X  Niunber  of  dete  polnte  for  frequency  varleble 
NW  -  1001; 

%  Define  number  of  polnte  per  2*pl/dw  for  time  domain 
NT  -  10; 

%  Temporal  llmlte  for  plotting  correlatlone 

teul  -  -6e-12; 

teu2  -  6e-12; 

teuhl  -  -6e-12; 

teuh2  ■  6e-12; 

%  Refractive  Index  of  5LT  et  40  C  (Bruner) 

ne.elt  -  Inline (['aqrt (4. 602483  +  (0.007294  +  3 . 483933e-8* (40  +  273.16)‘‘2) ./((2*pl*3e8./w*le6) .‘2  -  (0.185087  ♦  1. 607839e-8*(40  +  273. 15)‘2)‘2) 
'0. 073423. /((2*pl*3e6./w»le8). ‘2  -  0.199595'‘2)  +  0.001. /((2*pl*3e8. /w*le8)  . ''2  -  7.99724'*2)  -  0.02357*(2*pl*3e8./w*le6)  . ‘2)  *]  .  'w'); 

%  Refrectlve  Index  of  CLN  et  21  C(01eter  Jundt,  JOSAB  14,  3319)  0.4-5  urn 

ne.ln  -  lnllne([»eqrt(l  +  2.9804*  (2»pl*3e8. /w/le-8)  . '‘2. /U2*pl*3e8./w/le-6)  . ‘2  -  0.02047)  +  0.6961*(2*pl*3e8. /w/le-6)  .  "2.  / '  ... 

»((2*pl*3e8./w/le-«).“2  -0.0666)  ♦  8.9543*(2*pl*3e8./w/le-8) . "2./( (2*pl*3e8. /w/le-6) . "2  -416.08))’],  ’w*); 
no.ln  -  lnllne([»8qrt(l  +  2. 6734* (2*pl*3e6. /w/le-6) . ‘2. /((2*pl*3e6. /w/le-6) . “2  -  0.01764)  +  1. 2290* (2*pl*3e8. /w/le-6) . ‘2. / ’  ... 
’((2*pl*3e8./w/le-«).‘2  -0.05914)  +  12.614*(2*pl*3e8. /w/le-8) . “2./((2*pl*3e8. /w/le-6) . ‘2  -474.6))’],  ’w*); 

X  Refractive  Index  of  5%  HgO  doped  CLN  at  21  C  (Dieter  Jundt,  JOSAB  14,  3319)  0.4-5  urn 

ne_mgo61n  -  lnllne([’eqrt(l  ♦  2.4272*(2*pl*3e8./w/le-8) .'•2./((2*pl*3e8./w/le-6) .‘2  -  0.01478)  +  1.4617* (2*pl*3e8. /w/le-6) . ‘2. /*  ... 

’((2*pl*3e8. /w/le-6) .‘2  -  0.05612)  +  9. 6636* (2*pl*3e8. /w/le-6) . ‘2./ ((2*pl*3e8. /w/le-6) . ‘2  -  371.218))*],  ’w’); 
no_mgo51n  -  lnllne([*eqrt(l  +  2.2454*(2*pl*3e6. /w/le-6) .‘2./ ((2*pl*3e8. /w/le-6) .‘2  -  0.01242)  +  1.3006* (2*pl*3e6./ w/le-6) . ‘2. /*  ... 

’((2*pl*3e8. /w/le-6). ‘2  -0.05313)  +  6. 8972*(2*pl*3e8. /w/le-8) . "2./ ((2*pl*3e8./w/le-6) . "2  -331.33))*],  *w*); 

X  Refrectlve  Index  and  nonlinear  coefficient  for  generating  cryetal 
np  ■  ne.elt; 
ne  -  ne.elt; 
nl  ■  ne.elt; 

d  -  12.9e-12:  X  d33  Mg0:SLT 

X  Refrectlve  Index  and  nonlinear  coefficient  for  eummlng  cryetel 
npe  -  ne.ln; 
nee  ~  ne.ln; 
nle  -  ne.ln; 

de  -  26.0e-12;  X  d33  CLN 

X  Oeflne  pump  frequency 
wp  •  2*pl*c/632.10e-9; 

X  Oeflne  pxinp  power 
Pp  -  6.9; 

X  Define  elgnal  center  frequency  and  bandwidth 
weO  -  2*pl*c/1000e-9: 
dw  -  2*pl*3e 10*300; 

X  Define  the  QPH  order 
qpn.order  ■  1; 

X  Oeflne  average  duty  cycle 
dcycle  ■  0.5; 

X  Oeflne  domain  edge  poeltlon  nolee  (etandard  deviation,  fraction  of  Ldnolee) 

Xdnolee  -  l/eqrt(2)*0.375; 
dnolee  ■  0; 

X  Oeflne  domain  width  to  uee  for  domain  edge  poeltlon  nolee  fraction  (eet 
X  to  zero  to  use  actual  domain  width 
Ldnolee  ■  3.5e-6; 

XLdnolee  -  0; 

X  Define  mlnlmtim  etep  elze  of  mask  writer 
tmask  •  O.Oe-5; 

X  Set  to  1  to  uee  low  gain  assumption 
low.geln  ■  0; 

X  Oeflne  cryetel  length 
L  -  0.02; 
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X  S»t  to  1  for  cMrpod  cryetal 
chlrp«<igc  -  1; 

%  D«fln«  chirp  dlractlon  (radlaft  ■*  1  for  rad  at  z  ■■  0,  hlua  at  z  -  L) 
radlaft  -  0; 

X  Set  to  1  to  perfom  phaaa  coapaneatlon 
phasacoap  -  0; 

X  Sat  to  1  to  apodlza 
apodlza  ■  0; 

X  Sat  to  1  to  filter  the  output  epectrxin 
filtar_*eiii_epac  -  0; 

X  Define  relevant  fraquanclaa 
ve.hlgh  -  weO  ♦  dw/2; 
ire.low  -  weO  -  dw/2; 

we.hlgh.llnlt  -  2*pi*c/603.66a-9; 
we.low.llait  ■  2*pi*c/4500e-9; 

X  Sat  frequency  dooaln  for  3«dv  bandvldrh  cantered  at  vsO 
X  hut  additionally  hounded  hy  [vs.low.llnlt,  we_hlgh_llalt] 
wl  ■  wa_low  -  dw; 
w2  -  ve.high  ♦  dw; 

lf(we.hlgh  dw  >  ws.hlgh.lialt)  w2  ■  wa.hl^.lialt;  end; 
lf(we_low  -  dw  <  we_low_llait)  wl  -  w8_low_linit ;  end; 

X  Force  lower  bound  at  degeneracy  (algnal:  we  >  wp/2,  idler;  wp  -  we) 
if(wl  <  wp/2) 
wl  -  wp/2; 

and 

we  -  llnspace(wl,  w2,  NV) ; 
wl  -  wp  -  we; 

X  Define  frequency  interval  of  aein  epectrxm 

Ivlh  >  iBln(flnd(ve  >  we.lou))  :aax(f  lnd(we  <  we_hlgh)); 

xxxxxxxxxxxxxxxxxxxxxx 

X  Generating  Cryetal  X 
XXXXXXXXXXXXXXXXXXXXXX 

X  Coopute  k  vactora 
kp  ■  np(wp)*wp/c; 
ka  *  nsCwe) .*wa/c; 
ki  -  ni(wi).*wi/c; 

X  Compute  phase  mlematch 
Dk  -  kp  -  (ke  +  ki); 

if(chirpedgc  —  1) 

X  Compute  chirp  rate 
if(redleft  —  1) 

DkO  >  epllne(we,  Dk,  we.low); 

DkL  -  epllne(w8,  Dk,  we.high) ; 

alee 

DkO  -  eplina(w8,  Dk,  we.high); 

DkL  -  apline(w8,  Dk,  we.low); 

end 

zeta  -  -(DkL  -  DkO)/L; 

X  Compute  phase  matching  poeitions  for  linear  k-vector  chirp 
zpm  >  L  +  (DkL  -  Dk(iwlh))/zeta; 

alee 

DkO  -  epline(v8,  Dk,  weO); 
zeta  >  0; 
end 

X  Generate  vector  of  domain  lengthe 

dlep(aprintf ( ’\nGenaratlng  vector  of  domain  lengthe  (generating  cryetal)...*)); 

X  Define  chirping  function  (poling  period  as  a  fiinction  of  z) 
zO  -  lln8pace(0,  L,  1000); 

X  Calculate  group  valocitiae 
vgs  ■  der(we.  ke,  1); 
vgl  -  der(wi,  ki,  1); 
vgr  -  l./(l./vge  -  l./vgi); 

lf(chlrpedgc  »  i  kJt  phasecomp  —  1) 

X  Store  uncompeneated  phase  matching  poeition  vector 
zpm.orig  •  zpm; 

X  Compute  approximate  relative  (aignal  -  idler)  group  delay 
d]^l_ep  ■  (L  -  zpm.orig). /vgr(iwlh); 

X  Remove  third  order  dlepereion 

XP  ■  polyfit(w8(iwlh)  -  weO,  dphi.ep,  2); 

Xdphi.sp.pc  ■  dphi.ep  -  P(l)*(we(iwlh)  -  we0).“2; 
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%  R«(iOv*  fifth  ordsr  dispersion 

XP  ■  polyfit(vs(iwlh)  -  wsO,  dphi.sp,  4); 

%dphi_sp^pc  -  dphi.sp  -  P(i)*(ws(iwlh)  -  W80).''4; 

X  Force  linear  chirp 

Xdphi_8p_pc  •  (dphi_8p(end)  -  dphi_sp(i))/(v8_high  -  ws_low)*(ws(iwlh)  -  ws.low)  ♦  dphi_sp(i); 

X  Fores  zero  5th  and  higher  order  dispersion 
XP  “  polyfit(w8(iwlh)  -  wsO,  dphi.sp,  3); 

Xdphi.sp.pc  -  polyval(P,  ws(itflh)  -  wsO); 

X  Force  zero  3rd  and  higher  ordsr  dispersion 
XP  “  polyfit(w8(iwlh)  -  wsO,  dphi.sp,  i); 

X<^hi.8p.pc  -  polyvaKP,  wsCiwlh)  -  wsO); 

X  Fores  zsro  4th  and  higher  dispersion  and  fix  TOO/GOO  ratio 
P  -  polyfit(w8(iwlh)  -  wsO,  dphi.sp,  2); 

P(i)  -  i/2* (P(2)»3S9. 6/236. 8*le-i6);  X  2*P(i)  -  TOO,  P(2)  -  GOO,  ratio  used  here  is  for  fused  silica  at  966.2  nn 
dphi.sp.pc  -  polyvaKP,  wsCiwlh)  -  wsO); 

X  Ensure  that  pairs  phase  matched  at  ths  output  of  ths  crystal  havs  zsro 
X  relative  group  delay  (this  fixes  zpa  at  L  for  those  pairs) 
if(redlsft  —  i) 

dphi.sp.pc  -  dphi.sp.pc  -  spline (ws ( Iwlh) ,  dphi.sp.pc,  ws.hlgh) ; 

else 

dphi.sp.pc  -  dphi.sp.pc  -  splinsCtfsCiwlh) ,  dphi.sp.pc,  ws.low); 

snd 

X  Conputs  chirping  function 

zpm  ■  L  -  vgr(iwlh) . *dphi.8p.pc ; 

pp  *  2*pi./8pllns(zpa,  Dk(lwlh),  zO); 

elss 

pp  ■  2*pl./(DkO  -  zsta*zO) ; 

snd 

ifCtmask  0) 

X  Idsal  mask  writer 
Ld  -  []; 
z  -  0; 

whilsCz  <  L) 

dz  -  interpKzO,  pp,  z)*qpa.order*dcycl8; 
dcycle  ■  i  -  dcycle; 

Ld  -  [Ld  dz]; 

2  -  z  +  dz; 

snd 

elss 

X  Non-ideal  mask  writer 

z  -  0; 

za  >  0; 

zdr  -  []; 

whilsCzm  <  L) 

dz  -  interpKzO,  pp,  z)*qpn.order*dcycle; 
dcycls  ■  i  -  dcycle; 
z  -  z  +  dz; 

zm  -  zm  +  round(dz/tma8k)*tma8k; 
if(z  -  zm  >  tmask/2) 
zm  >  zn  +  tmask; 

else 

if(z  -  zn  <  -tma8k/2) 
zm  -  za  -  tmask; 

snd 

snd 

zdr  ■  [zdr  ai] ; 

snd 

Ld  -  diff([0  zdr]); 

snd 

X  Add  domain  noise 

X  Fora  a  vector  of  thu  positions  of  ths  right  edges  of  each  domain 
zdr  -  cuasua  (Ld) ; 

X  Perturb  thess  edges 
ifCLdnoiss  —  0) 

for  ii  -  i: length (Ld) 

zdr(ii)  >  zdr(il)  dnoi8e*Ld(il)*sqrt (-2*log(rand))«co8 (2*pi*rand) ; 

end 

slss 

for  il  -  i; length (Ld) 

zdx(ii)  >  zdr(il)  *  dnol8e*Ldnoise*sqrt(-2«log(rand))*co8(2*pl*rand) ; 

end 

end 

X  Reconstruct  domain  length  vector 
Ld  -  diff([0  zdr]); 

X  Fix  crystal  length 
L  •  8um(Ld): 

X  Perform  apodizing 

apodaask  •  onss(i,  Isngth(Ld)); 

lf(chirpedgc  —  i  tt  apodize  --  1) 
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%  D«fiae  spodlzstlon  function  In  polynoolsl  form 

%  Tbs  apodlzatlon  function  should  be  tbs  probability  that  s  rsvsrssl  is 
%  KOT  rsnovsd. 

XPapod  -  [-5. 06930464 14 99887 s-45  2. 84946298 ii0ii68s -29  -6. ii6363656404832s-i4  29.64933346526240]; 

%fspod  ■  epllnsCtpa,  polyvsKPspod,  vs),  zdr) ; 

%  Assums  sll  frequency  components  confocslly  focused 
Ap  ■  2*pi*c/wp*L/(4*np(wp)) ; 

As  -  2*pi*c./v8*L./(4*ns(ws)); 

Al  -  2*pi*c./wl*L./(4*ni(wi)); 

etap  ■  ets0/np(wp); 
etas  -  stsO. /nsCvs) ; 
stsl  -  stsO./ni(vl); 

%  Computs  gsln  constant 

kappa  >  eps0*d*sqrt(8*hb*wp*v8.*vl«'etsp.«etas.  *etal*Ap.*As.*Ai  ... 

./(Ap*As  ♦  Ap*Ai  ♦  Aa.*Ai).‘2); 

fspod  "  kappa(lvlb(8nd) ) . /kapps(ivlh) ; 

X  Computs  spodlzstlon  function  (probability  vs.  z  thst  s  rsvsrssl  le 
\  NOT  removed) 

fspod  "  splineCzpm,  fspod,  zdr); 
wlnf  -  3; 

win  >  (1  -  cos(2*pl* (l;round(l8ngth(fspod)/wizif ))/(length(f 8pod)/winf )) )/2; 

vln  •  [vln(l:round(length(vin)/2))  onesCl,  lengthCfspod)  -  length(win))  win(round(length(wln)/2)>l:end)] ; 
Xfapod  ~  fspod. *wln; 

X  Fix  erroneous  downturn  due  to  polynomlsl  end  behavior 
XfspodCfind (fspod  —  Dax(f8pod)) ;end)  ■  i; 

K  Compute  target  fitnction  to  compare  against  when  integrating 
%  normalized  nonlinearity 
target  >  cunsimp(zdr ,  fspod); 

%  Computs  period  widths 
if  (mod(l8ngth(Ld)  .  2)  —  0) 

Lp  -  Ld(i:2;end)  +  Ld(2:2:end); 

else 

Lp  -  Ld(l;2: end-1)  *  Ld(2:2:end); 
end 

'f,  Computs  period  mask  (0  for  no  reversal,  1  for  reversal) 
prmask  -  onesd,  floor(length(Ld)/2)) ; 

for  il  -  1; length (Lp) 

if(8Us(Lp(l:ll) .*prmask(l:il))  >  targ8t(ll*2  -  i)) 
pmaskdi)  >  0; 

snd 

snd 

8podma8k(2:2;end)  ■  2*prm8sk  -  i; 
end 

mxnmrmrmw. 

%  Summing  Crystal  % 

mxnmmmxxxx 

%  Compute  k  vectors 
kpe  ■  np8(wp)*wp/c; 
kss  -  ns8(vs) .*us/c; 
kls  -  nl8(wl).*wl/c; 

X  Compute  phase  mismstch 
Dks  >  kps  -  (kss  *  kle); 

X  Compute  phase  mismstch  at  center  frequency 
DksO  -  epllns(v8,  Dks,  vsO) ; 

X  Modify  crystal  length  as  necessary  to  achieve  s  doss  match  correlation 
Ls  -  i8-3; 

X  Gensrsts  vector  of  domain  lengths 

di8p(8printf ('Generating  vector  of  domain  lengths  (summing  crystal)...')); 

X  Define  chirping  function  (poling  period  as  s  function  of  z) 

zO  •  lln8pacs(0,  Ls,  1000); 

pp  -  onssd,  lsngth(z0))*2*pi. /DkeO; 

lf(tmask  “  0) 

X  Ideal  mask  writer 
Lds  -  []; 
z  -  0; 

vhlle(z  <  Le) 

dz  -  Int8rpl(z0,  pp,  z)*qpm_ord8r*dcycle; 
dcycls  ■  1  -  dcycle; 

Lds  -  [Lds  dz] ; 
z  -  z  +  dz; 

end 

slss 
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%  Non-id«al  mask  uritsr 
z  -  0; 
za  -  0; 
zdrs  -  [] ; 
wliilsCza  <  Ls) 

dz  •  int«zpi(z0,  pp,  z)*qpo_order*dcycle; 
dcycls  ■  1  -  dcycla; 
z  -  z  ♦  dz; 

zm  ^  zm  *  round (dz/toask)*taa8k; 
if(z  -  zm  >  taask/2) 
zm  •  zm  *  tmask; 

els* 

if(z  -  za  <  -tmask/2) 
za  •  za  -  taask; 

•nd 

end 

zdrs  -  [zdre  zm] ; 

end 

Lds  -  diff([0  zdrs]); 

end 

X  Add  dcnain  nois* 

X  Fora  e  vector  of  the  position*  of  the  right  edge*  of  each  domain 
zdrs  -  cumsuaCLde) ; 

X  Perturb  these  edges 
if(Ldnoiss  —  0) 

for  11  •  1; length (Lda) 

zdrs(il)  *  zdre(ii)  +  dnoiee*Lds(ii)*sqrt(-2*log(rand))*co8(2*pi*rand) ; 

end 

ele* 

for  11  *  i:length(Lds) 

zdra(ii)  •  zdre (11)  +  dnoi8e*Ldnoise*sqrt(-2*log(rand) )*C08(2*pi*rand) ; 

end 

end 

X  Reconstruct  domain  length  vector 
Lds  ■  diff([0  zdre]); 

X  Fix  crystal  length 
La  ~  aum(Lds); 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X  Compute  ABCD  Matrix;  Generating  Crystal  X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X  Compute  aqrt  pump  photons/s 
bp  -  eqrt(Pp/(hb*wp)) ; 

X  Assume  all  frequency  components  confocally  focused 
Ap  «  2*pl*c/wp*L/{4*np(wp)) ; 

As  -  2*pi*c. /ws*L./(4*ns(ws)) ; 

Ai  ■  2*pi*c. /wi*L./(4*ni(wi)) ; 

step  «  eta0/np(vp); 
etas  -  ets0./ns(us) ; 
etai  >  etaO. /ni(vi) ; 

X  Compute  gain  constant 

kappa  *  apa0*d*8qrt(8*hb*wp*ws . *ui*etap. *etas. *atai*Ap. *Aa. *Ai  ... 

./(Ap*Aa  +  Ap«Ai  +  Aa.*Al)  . '‘2)  ; 

X  Compute  overlap  factors  separately  for  reference 
ge  «  2«Ap.«Ai./(Ap.*Aa  Ap.*Ai  +  As.*Ai); 
gi  -  2*Ap.*As,/(Ap.*As  ♦  Ap.*Ai  +  As.*Ai); 
gp  “  2*Aa.*Ai./(Ap.  «Aa  Ap.*Ai  +  As.*Ai); 

dlsp(sprintf ('Computing  ABCD  matrix  (gensrsting  crystal)...*)); 

Al  -  zeroed,  langth(we)); 

Bi  -  zerosd,  length(we)); 

Ci  -  zerosd,  langth(we)); 

Di  -  zeroed,  length(we)); 
zdl  -  cumsum([0  Ld(l;snd-1)] ) ; 
egnd  >  2*mod(l: length (Ld) ,  2)  -  i; 
sgnd  -  sgnd.*spodmask; 
lf(low_galn  ■«  1) 

X  Low  gain  assumption:  simply  integrate  signal  and  idler  propagation 
X  equations  independently 
for  il  -  i:length(vs) 

Mi2  -  sum(sgnd*ks]^s(il)*bp/Dk(il) . *(ezp(i*Dk(il)*(zdl  ♦  Ld))  -  *xp(i*Dk(il)*zdl>) )  ; 
Aidl)  -  exp(i*ks(ii)*L); 

Bidi)  -  M12**xp(i*ka(ll)*L)  ; 

Cidi)  -  conJ(Mi2)*exp(-i*ki(ii)*L); 

Di(ii)  -  *xp(-i*ki(ii)*L); 

and 

else 

X  High  gain:  integrate  coupled  equations  for  each  slab  iteratively 
X  Note  that  only  the  sign  of  keppe  changes  domaln-to-domein;  e  ■  const 
s  •  sqrt(bp“2*ksppa. “2  -  Dk.“2/4); 

Ali  -  sxp(i*Dk*Ld(l)/2)  .•(cosh(8*Ldd))  -  i*Dk./(2*8) .  •8inh(8*Ld(l))) ; 

Bii  ■  i*bp*kepp*./e.**xp(i*Dk*(zdl(i)  +  Ld(l)/2)) .*ainh(8*Ld(l)) ; 
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for  11  ■  2: length (Ld) 

A21  -  exp(l*Dk*Ld(ll)/2>.*(co«h(8*Ld(ll))  -  l*Dk./(2*8)  .  *8lnh(e*U(ll)))  ; 

B21  -  l*bp*8gnd(ll)*k8pp*./8.*8Xp(l*Dk*(zdl(ll>  ♦  Ld(ll>/2)) .♦8lnh(8*Ld(li) ) ; 

Alin  -  A11.*A21  ♦  conj(Bll) .*621; 

Blln  -  B11.*A21  ♦  conj(All>.*B21; 

All  -  Alin; 

Bll  •  Blln; 

end 

Al  •  All  .*exp(l*]iC8*L)  ; 

Bi  •  Bll. *exp<l*k8*L) ; 

Cl  •  conj (Bll) . •exp(-l*kl*L) ; 

Di  -  conj (All) .•8xp(-l*kl*L) ; 

snd 

X  Apply  a  band  pans  filter  to  tbs  asln  spectrum 
If (filter jBsln. spec  “  1) 

HF  ■  onesd,  l8ngtb(w8))  .•(ws  >  ws.low  -  dw/4  A  we  <  ws.blgh  ♦  dw/4) ; 

Al  -  A1.*HF; 

Bl  -  B1.*HF: 

Cl  -  Cl . ♦conj (HF) ; 

Di  -  Di.*conJ(HF); 

snd 

mmxxxxmxxxxxxrmxxnmxxnmmmx’xx 

X  Coaputs  Transfer  Function;  Sunning  Crystal  X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X  Assuds  all  frequency  components  confocally  focused 
Ape  ■  2»pi*c/wp*L8/(4*np8(wp) ) ; 

Ass  *  2*pl*c./w8*Ls./(4*nse(ve)) ; 

Ais  ■  2*pl*c./wi*Ls./(4*nl8(wl)) ; 

etaps  >  ets0/np8(wp) ; 
stass  *  8t80./n8s(w8) ; 
etais  -  etsO. /nls(wl) ; 

X  Compute  gsln  constant 

kj^>pas  ~  epsO^ds^sqrt (8*hb*vp* vs wiretaps. ♦stsss . *8tai8*Ap8 Ass.  * Ais  ... 

./(Aps^Ass  Aps*Als  *  Ass.^Als)  .*2) : 

X  Compute  overlap  factors  separately  for  refsrencs 
gss  -  2*Ape.«Ais./(Aps.*Ass  *  Ape.*Als  *  Ass.*Als); 
gis  -  2*Ape.*A88./(Aps.«Ass  *  Ape.* Ais  *  Ass.*Als); 
gpe  >  2*A8S.*Als./(Ap8.*Ae8  *  Ape.* Ale  *  Ass.*Als); 

dlsp(sprlntf ('Computing  transfer  fxmctlon  (summing  crystal)...')); 

ESP  *  zeroed,  l8ngtb(u8)); 
zdls  -  cum8um([0  Lds(l; snd-i)]) ; 
sgnds  -  2*mod(l : length ( Lds) ,  2)  -  1; 
for  11  -  l;length(ws) 

M12  -  sum(sgnd8*kappss(li)/Dks(il) .*(exp(-i*Dks(ll)*zdls)  -  sxp(-l*I^8(ll)* (zdls  +  Lds)))) 
ESP(ii)  -  Mi2*exp(i*kpe*L); 

end 

XXXXXXXXXXXXXXXXX 
X  PSD  and  Pbass  X 
XXXXXXXXXXXXXXXXX 

X  Compute  signal  power  spectral  density 
Ssl  -  i/(2*pi)*abe(Bl).-2; 

X  Calculate  paired  gsnsrstion  rets 
R  >  simp(v8,  Ssl); 

X  Plot  signal  powsr  spectral  density 
f  ig\ire: 

plot (vs,  Ssl,  'k-'); 

t It le(8prlntf ('Signal  power  spactrsl  density  (R  ■  X.3e  (8"{-l}))',  R)); 
xlabsK 'Xomsgs  (rsd/s)'); 
ylabel('S_s(\omega)  (8“<-l>/rad) ') ; 

X  Compute  blphoton  phase  (signal  *  idler),  relative  delay,  and  GI^ 
phi  -  unurap(angls(Al.*conJ(Cl)),  pi); 
dphl  -  der(phl,  vs,  1); 
d2pbl  -  der(phl,  ws,  2); 

X  Plot  phase,  delay,  and  CDD 

figure; 

subplot (311) ; 

plot ( vs ( ivlh) ,  phi ( Iwlh) ,  ' k- ' ) : 
axis  tight; 

tit Is ('Biphoton  phase*) ; 
xlsbslCXomsga  (rsd/s)'); 
ylabeK 'Phase  (rad)'); 
subplot (312) : 

plot (vs (iwlh) ,  dphl (Ivlh),  'k-'); 
axis  tight; 

title ('Signal/Idler  relative  group  delay'); 
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xl«b«l(’\o<B«ga  (rad/s)*); 
ylab«l(*Group  daisy  (s)’): 
subplot (313) : 

plot(W8(iwlh),  d2phi(iwlh),  'k-'); 
axis  tight; 

tltloC 'Blphoton  GDD’); 
xlabsK  ’  \oaaga  (rad/s)  * )  ; 
ylabaK’GDD  (8*2/rad)'); 

%  Calculats  avaraga  GDD  and  TDD  based  on  nth  ordsr  fit  to  phass  In  significant 
%  fraquancy  rsgion 

P  -  polyf it(w8(iwlh)  -  wsO,  phl(lwlh),  3); 

GDD  -  2*P(2); 

TOD  -  6*P(i): 

*^D6D  -  i20*P(i); 

%  Show  fit  in  phass  plot 
%subplot(3ii) ; 

Xhold  on; 

%plot(ws(iwlh)  ,  polyvaKP,  wsClvlh)  •  wsD)  ,  'k — ’); 

X  Spaclfy  Optical  Systan  and  Coapraaaor  X 

xxxxxmxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X  Signal  optical  system  transfer  function 
HI  -  1; 

X  Idlsr  optical  system  transfsr  function 
G1  -  1; 

X  Compressor  transfer  function  (signal  path) 

XHC  -  l; 

XHC  -  axp(-i*phi); 

X  Pentolsssr  costing 

phic  -  -4.64Da-42*ws."3  +  2.662a-26*ws. ‘2  -  D. 00226*sin(3. 3333s- 13*ws) ; 
phlc  •  phic  -  polyv al (polyf it (ws.  phic,  i) ,  ws);  X  Raaova  time  daisy 

phic  -  phic* 13;  X  Number  of  bounces 

HC  •  axp(l*phic) ; 

X  Compressor  transfer  function  (idlsr  path) 

GC  -  i; 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X  Uncompressed  Correlations  X 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

di8p(sprlntf (’Computing  uncomprsssed  corrslations ...’)) ; 

X  Spaclfy  ABCD  coefficients  incident  on  the  summing  crystal  (i.a.  after  ths 
X  optics!  system) 

A2  -  A1.*H1; 

B2  -  Bl.*Hi; 

C2  -  Cl. *conj (Gl) ; 

D2  -  Dl.*conj(Gl); 

X  Use  Fourier  transform  to  compute  vavef unction 
Ws  -  FTCraata(w8(2)  -  wad),  2*pl/dw/NT) ; 
tsu  -  FTDcmain(W8) : 

X  Compute  the  uavefunction  without  summing 

p8i2  *  i/(2*pi)*exp(i*(up  -  w8D)*t8u) . *IFTCP(interpi(w8  -  wsD,  A2.*conJ (C2) ,  Ws,  ’nsarsst’,  D) ,  Ws) ; 
fhl  -  ehs(p8i2)  .*'2; 

X  Compute  the  uncompressed  SFG  rats 

Rmq  -  l/(4*pi“2) .*8b8(IFTCP(interpl(ws  -  weD,  A2.*conJ (C2) .*ESP,  Wa,  ’nearest',  0),  W8)).“2; 

X  Plot  wavsfunction  amplitude  squared 
figure; 

plot(tau,  fhl,  'k-’); 

V  -  aria;  axi8([t8ui  t8u2  V(3)  V(4)]); 

tltle(’Biphoton  intensity  (uncompressed)  l<a_8(t)  a.i(t  -  \teu)>|"2’); 
zlebel(’\teu  (a)’); 

X  Plot  quantum  sum  photon  rats  for  the  non-chirped  crystal 
figure; 

plot(tau,  Rmq,  'k-',  tau,  fhi/max(fhl)*maz(Rmq) ,  ’k — ’); 

V  -  axis;  axl8([t8ui  tau2  V(3)  V(4)]): 

title(8printf( 'Non-chirped  quantum  aum  photon  rate  (uncompressed,  L_8  ■  X.3f  (on))',  Ls/le-S)); 
zlabel(’\tau  (a)'); 
ylabel('R.{mq}  (8“{-l»’); 

legend(’R_'(mq}* ,  '\propto|<e_8(t)  e.l(t  “  \tau)>|“2',  'Location',  'Best'); 

xxxxxxxxxxxxxxxxxxxxxxxxxxx 

X  Compressed  Corrslations  X 
XXXXXXXXXXXXXXXXXXXXXXXXXXX 

dlep(sprlntf( 'Computing  compressed  correlations...’)); 

X  Specify  ABCD  coefficients  incident  on  the  suimnlng  crystal  (i.s.  after  the 
X  optics!  system) 
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A2  - 
B2  - 

C2  -  Cl.*coiij(Gl)  .•conjCGC); 

D2  -  Dl.*conj(Gl).*conJ(GC); 

%  Computs  ths  wavsfimction  without  sunning 

p8i2h  •  l/(2*pi)*«xp(i*(wp  -  w8D)*tau)  .*IFTCP(int8rpl(w8  -  W8D,  A2.  *conJ(C2) ,  Vta,  *n«»r«8t’,  0),  We); 
fbih  «  abs(p8l2h)  .*'2; 

%  Compute  ths  uncomprsssed  SFG  rate 

Roqh  -  l/(4*pi"2>  .•sbsCIFTCPCinterpKws  -  wsO,  A2.  *conj  (C2)  .  *ESP,  Ws,  'nearest',  D)  ,  W8))."2; 

%  Plot  wav8f\mction  amplitude  squared 
fig\irs; 

plotCtau,  Ibih,  'k-’>; 

V  -  axis;  axisCCtauhl  t8uh2  V(3)  V(4)]); 

titlsCBiphoton  intsnsity  (comprsesed)  |<s_8(t)  a_i(t  -  \t8u)>l"2'>; 
xlsbel(’\tau  (s)’); 

%  Plot  quantum  sum  photon  rsts  for  the  non^chirped  crystsl 
figure ; 

plotCtsu,  Rmqh,  'k-*,  tsu,  fbib/max(fbib)*max(Rmqb) ,  'k — '); 

V  -  axis;  axisdtsuhl  tauh2  V(3)  V(4)]); 

titlsCsprintf ('Non-cbirped  quantum  sum  photon  rsts  (comprsssed,  L.s  -  %.3f  (mm))',  L8/le-3)); 
xlsbeK’Xtau  (s)’); 
y label (’R_{mq>  (8“-(-l»'); 

leg8nd( ’R_{mq}’ ,  '\propto|<8_8(t)  8_i(t  -  \t8u)>|“2',  'Location',  'Best'); 

%  Combine  plots  onto  one  figure 

figure; 

subplot (221); 

plot (2*pi*c./w8» 189,  Ssl,  'k-'); 

V  -  axis;  axi8( [2*pi*c/wa(and) *109  2*pi*c/wa(l)*l89  V(3)  V(4)]); 

titl8(-{' Signal  power  spectral  density*;  8printf('(R  ■  X.3e  (8‘'{'1}))',  R)}) ; 

xlsbeK 'Mambds  (nn)'); 

ylsbeK 'S_s(\lambd8)  (8“-(-l}/r8d)  ') ; 

subplot (222) ; 

plot (2*pi*c./w8(iwlh) *189,  phi(iwlb),  'k-'); 
axis  tight; 

titlsCBiphoton  phasa  spectrun'); 
xlsbsK'Mambda  (nm)'); 
ylabsK'Pbass  (rad)'); 
eubplot (223) ; 

plot(tau,  Rmq,  'k-',  tsu,  fbi/max(fbi)*max(Rmq) ,  'k — '); 

V  -  axis;  axi8([t8ul  t8u2  V(3)  V(4)]); 

title ({'Quantum  SFG  rsts';  sprintf(' (uncomprsssed,  L_s  ■  (on))',  Ls/ls-S)}); 

xlsbeK '\tsu  (s)*): 
ylabel('R_{mq>  (8“{-l})'); 

leg  «  legend ( 'R.'Cmq}' ,  ’\propto| <8_8(t)  8_i(t  -  \tau)>|“2',  'Location',  'Best'); 
setdeg,  'FontSiza',  ID); 
subplot (224) ; 

plot(tsu,  Rmqh,  'k-' ,  tsu,  fbih/max  ( fbih)  *m  ax  (Rmqh)  ,  'k — '); 

V  «  axis;  axi8([t8uhl  t8uh2  V(3)  V(4)]); 

tit  Is ({'Quantum  SFG  rata';  sprintf (' (compressed,  L.s  ■  X>3f  (mn))',  Ls/l8-3)}); 
xlsbeK *\tau  (s)'); 
yl8b8K'R_{mq>  (a“{-l})'); 

leg  «  leg8nd('R_{mq}' ,  '\propto|<8_8(t)  a_i(t  -  \t8u)>l“2',  'Location',  'Bast'); 
sstdsg,  'FontSize',  ID); 

%  Report 

disp(8printf ('VnPairsd  rsts:  %.3s  (s''-!)',  R)); 

di8p(8printf('GDD  (based  on  order  (X.Of)  fit):  X.Ss  (f8“2)',  l8ngth(P)  -  1,  GDD/le-30)); 
di8p(sprintf ('TDD  (based  on  order  (X.Df)  fit):  X.3e  (fs"3)',  langth(P)  -  1,  TOD/la-46)); 
disp(eprintf ( 'Non-chirped  summing  crystal  length:  X.3f  (an)',  Ls/ls-3)); 

if(cbirpedgc  —  1) 

di8p(8printf ('Poling  period  (generating  crystal):  Rad:  X.4f  (urn).  Blue:  X.4f  (xn) ' ,  ... 

2*pi/splins(ws,  1%,  ws.low)*qpm_ordsr/ls-6,  2*pi/8plina(w8,  Dk,  w8_high)*qpQ_ordar/la-6) ) ; 

slss 

dl8p(8printf ('Poling  psriod  (generating  crystsl):  X.4f  (m)*,  2*pi. /DkO*qpQ_ord8r/l8-6)) ; 

snd 

disp(eprintf (’Poling  period  (siuning  crystsl):  X.3f  (um)  ' ,  2*pi. /Dk8D*qpn_ord8r/ls-6)) ; 

X  Display  poling  periods  for  SLT  st  25  C 
TsxpandSLT  -  1  +  1.68-5* (40  -  25)  ♦  78-9* (40  -  25) “2; 
if(chirpedgc  ■■  1) 

disp(8printf ('Poling  psriod  (generating  crystsl)  st  25  C  (SLT):  Red:  X.4f  (in).  Blue:  X.4f  (um) ’ ,  ... 

2*pi/8pline(w8,  Dk,  w8_low)/TexpandSLT*qpm.ordsr/ls-6,  2*pi/8plin8(ws,  Dk,  ws_hi^)/T8xpandSLT*qpffl_ord8r/l8-6)) ; 

else 

di8p(8printf ('Poling  period  (generating  crystal)  st  26  C  (SLT):  X.4f  (um) ' ,  2*pi./DkO/TaxpandSLT*qpo_ord8r/l8-6)) ; 

snd 

X  Display  poling  periods  for  CLK  st  25  C 
TsxpandCLK  -  1  14.8e-6; 

disp(8printf( 'Poling  period  (sunning  crystal)  at  25  C  (CLN) :  X.4f  (\n) ' ,  2*pi ./DkO/TaxpandCLN*qpm.order/la-6) ) ; 
slspsedtims  •  toe; 

di8p(8prlntf ('VnElspsed  time:  X-lf  (s)',  elapsedtine)); 
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Chapter  4 

Theory  of  EIT  based  Paired 
Photon  Generation 


Correlated  and  entangled  photon  pairs  are  widely  used  in  quantum  communication^ 
quantum  cryptography  [31]  and  quantum  imaging  [32,  33].  These  photon  pairs  are 
usually  produced  via  spontaneous  parametric  down  conversion  in  nonlinear  crystals. 
A  few  years  ago  a  new  approach  to  the  generation  of  paired  photons  has  been  exper¬ 
imentally  demonstrated  by  two  groups,  both  use  electromagnetically  induced  trans¬ 
parency  (EIT)  to  generated  paired  photons  in  an  otherwise  opaque  atomic  medium. 
Working  with  hot  atoms,  Lukin  and  colleagues  have  demonstrated  correlation  be¬ 
tween  generated  pulses  of  light,  as  well  as  storage  and  delayed  extraction  [34].  Work¬ 
ing  with  a  MOT,  Kimble  and  colleagues  have  shown  the  generation  of  nonclassical 
photon  pairs  with  a  programmable  delay  [35].  Recently,  the  Harris  research  group  at 
Stanford  has  demonstrated  generation  and  rudimentary  waveform  control  of  narrow 
band  biphotons  [36].  More  recently,  Kolchin  and  colleagues  have  shown  paired  pho¬ 
ton  generation  with  a  single  pump  beam  in  a  right  angle  geometry  and  paired  photon 
generation  in  the  ensemble  of  cold  two- level  atoms  [37].  A  long  coherence  length  and 
a  controllable  bandwidth  of  the  generated  paired  photons  are  the  advantages  of  the 
new  approach,  which  might  be  useful  for  such  applications  as  long  distance  quantum 
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communication  [38]  as  well  as  biphoton  waveform  control  and  shaping  ^ 


4.1  Theory  of  paired  photon  generation 


This  chapter  describes  the  theory  of  correlated  paired  photon  generation  in  a  collec¬ 
tion  of  double  A- type  atoms  [39].  Using  the  Heisenberg-Langevin  method  we  evaluate 
and  analyze  spectral  characteristics  of  the  generated  Stokes  and  anti-Stokes  photons 
and  their  time-correlation  properties.  Of  importance,  we  predict  the  regime  when 
the  correlation  time  of  generated  photons  is  determined  by  the  group  delay  caused 
by  EIT.  This  theoretical  prediction  as  well  as  some  theoretical  aspects  of  this  work 
has  been  outlined  in  Ref.  [36].  The  complete  theoretical  treatment  of  paired  photon 
generation  with  and  without  the  approximation,  that  the  atomic  population  remains 
in  the  ground  state,  has  been  given  in  Ref.  [39]. 


Figure  4.1:  Energy  level  diagram  and  schematic  for  spontaneous  backward- wave 
paired  photon  generation  in  an  atomic  cloud  formed  by  double- A  type  atoms.  In 
the  presence  of  the  pump  and  coupling  lasers  phase-matched,  counter-propagating 
Stokes  and  anti-Stokes  photons  are  generated  into  opposite  directions. 

The  schematic  of  the  process  considered  here  is  shown  in  Fig.  5.7.  In  the  presence 
two  cw  beams  termed  as  the  pump  and  coupling  lasers  with  frequencies  cUp  and 

^In  a  right-angle  geometry,  where  Stokes  and  anti-Stokes  photons  are  generated  and  collected 
at  right  angles  from  the  direction  of  the  pump-coupling  axis,  applying  the  absorption  mask  on  the 
pump  beam  allows  us  to  create  a  bi photon  with  a  precribed  waveform. 
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paired  spontaneous  photons  termed  as  Stokes  and  anti-Stokes  are  generated  in  the 
atomic  cloud  and  propagate  in  opposite  directions  along  the  2  axis.  In  order  to  keep 
the  parametric  gain  small  we  choose  the  pump  beam  to  be  weak  and  detuned  from  the 
resonance  transition  |1)  ^  |4).  The  intense  coupling  beam  is  tuned  to  resonance  with 
the  1 2)  ^  1 3)  transition  to  enhance  the  atom- field  interaction  and  provide  EIT  for  the 
generated  anti-Stokes  photon.  Under  such  conditions  we  expect  a  small  fraction  of  the 
atomic  population  to  be  in  the  excited  states  while  most  of  it  remains  in  the  ground 
state  |1).  The  frequencies  of  the  generated  photons  obey  the  energy  conservation 
+  ^as  =  +  ^c-  Ill  presence  of  EIT  the  anti-Stokes  photon  escapes  out  of  the 

atomic  cloud  with  very  slow  group  velocity. 

We  note  the  connection  to  earlier  work:  Two-photon  entanglement  in  type-II 
SPDC  has  been  analyzed  [40].  The  possibility  of  quantum  correlated  and  squeezed 
fields  in  backward  wave  EIT  system  has  been  predicted  [41]  and  large  parametric 
gain  and  oscillations  have  been  observed  [42].  Control  of  single  photons  has  been 
discussed  [43].  We  also  note  the  early  studies  on  double- A  atomic  systems  [44]. 

In  this  chapter  we  review  and  discuss  in  detail  the  system  dynamics  putting  em¬ 
phasis  on  the  infiuence  of  the  EIT  window  on  paired  photon  generation  bandwidth. 
We  show  that  at  low  parametric  gain  the  atomic  system  can  operate  in  two  different 
regimes.  In  the  first  regime,  where  the  group  delay  is  small,  the  intensity  correlation 
function  shows  Rabi  oscillations.  In  the  second  regime,  where  the  optical  depth  of  the 
atomic  sample  and  the  group  delay  are  large,  phase-matching  becomes  the  dominant 
process  that  controls  the  shape  of  the  intensity  correlation  function.  We  examine 
the  conditions  required  for  the  system  to  operate  in  the  oscillatory  and  group  delay 
regimes  and  discuss  corresponding  Stokes  and  anti-Stokes  spectral  generation  rates. 
In  particularly,  we  predict  that  the  group  delay  regime  requires  the  optical  depth 
much  higher  that  10.  This  point  has  been  missed  in  our  earlier  publication  [36]. 

We  also  extend  our  theoretical  treatment  of  paired  photon  generation  to  go  beyond 
the  ground  state  approximation.  This  allows  us  to  properly  include  and  analyze  the 
effect  of  Langevin  noise  fluctuation  on  the  atomic  system  and  solve  the  problem  of  its 
return  to  the  ground  state  after  the  emission  of  Stokes  and  anti-Stokes  photon  pairs. 
We  also  address  the  important  questions:  I)  “What  are  the  effects  of  the  optical 
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thickness  of  the  atomic  sample  on  paired  and  single  photon  generation?”  2)  “Does 
every  Stokes  photon  have  its  paired  anti-Stokes  photon?” 


We  introduce  and  derive  the  biphoton  wavefunction  taking  into  account  Langevin 
noise  terms.  We  also  derive  and  discuss  the  relation  of  the  Stokes  photon  count  rate 


to  the  classical  gain  coefficient  in  the  Stokes  channel. 

4.2  Heisenberg-Langevin  description  of  paired  pho¬ 
ton  generation 

We  consider  a  collection  of  identical  double  A-type  atoms  uniformly  distributed  within 
a  pencil-shaped  volume  with  cross  section  S  and  length  L.  We  assume  that  the  atomic 
sample  is  optically  thin  in  the  transverse  direction,  so  that  there  is  no  radiation 
trapping  effect  in  this  direction.  No  restrictions  are  imposed  on  the  optical  thickness 
of  the  atomic  sample  in  the  ^  direction.  We  also  assume  that  the  pump  and  the 
coupling  beams  counter-propagate  undepleted  through  the  atomic  medium.  Under 
these  assumptions  we  consider  propagation  of  a  single  transverse  spatial  mode  of 
radiation  along  the  z  axis.  The  pump  and  coupling  laser  beams  are  treated  as  classical 
quantities  and  their  interaction  with  the  medium  is  described  semi-classically.  In  order 
to  allow  for  the  spontaneous  initiation  of  the  parametric  fluorescence  process,  the 
generated  weak  Stokes  and  anti-Stokes  fields  are  described  by  quantum-mechanical 
operators,  in  slowly  varying  envelope  approximation: 


(4.1) 


where  subscript  j  denotes  either  Stokes  or  anti-Stokes  photon,  lJs  =  —  U2  +  Aa;i4, 

Y  =  L  X  S  is  the  interaction  volume. 

Adopting  the  notations  of  Lukin  and  Fleischhauer  [45,  46],  in  the  rotating  wave 
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approximation  we  write  the  interaction  Hamiltonian  in  continuous  form  as 


l/  = 


t)  +  gz\aas{z,  t)au{z,  t) 


+  y^32(2,  t)  +  ^^41(2,  t)  +  ^42^5(2, 0^42(2,  t)  +  h.C. 


(4.2) 


w'here  d'jk{z,  t)  are  the  collective  slowly  varying  atomic  operators,  defined  in  the  Ap¬ 
pendix  4.7.1,  N  is  the  total  number  of  atoms  in  the  atomic  ensemble,  Awh  is  the 
detuning  of  the  pump  laser  from  the  |1)  — »•  |4)  transition,  Qc  = 
are  pump  and  coupling  laser  Rabi  frequencies  with  Ep,Ec&s,  the  complex  amplitudes 
of  the  electric  fields,  Qs  =  and  gas  =  )f°*  are  the  coupling  constants  with  pjk 

as  the  dipole  moment  for  the  |j)  — >•  \k)  transition  and  Sj  =  as  the  electric  field 

of  a  single  photon. 

The  propagation  of  the  Stokes  and  anti-Stokes  fields  and  their  interaction  with 
the  atoms  are  described  by  the  set  of  Maxwell  and  Heisenberg-Langevin  equations. 
The  Heisenberg-Langevin  equations  are  responsible  for  the  atomic  evolution: 


=  +  rfk  +  Fjk, 


(4.3) 


where  ^jk  are  the  dephasing  rates,  are  the  spontaneous  emission  rates,  Fjk{z^t) 
are  the  collective  atomic  5-correlated  Langevin  noise  operators.  The  full  set  of  the 
Heisenberg-Langevin  equations  is  shown  explicitly  in  Sec.  4.7.2 

The  fluctuations  of  5-correlated  collective  Langevin  noise  operators  Fjk{z^t)  is 
given  by 

{Fjk{z,t)Fj,k'{z',t'))  =  ^'Djkj'k'{z,t)5{t  -  t')5{z  -  2'),  (4.4) 

where  Djkj'k'  is  a  Langevin  diffusion  coefficient.  The  derivation  of  Eq.(4.4)  and 
relevant  diffusion  coefficients  is  shown  in  Sec.  4.7.4. 
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The  evolution  of  the  annihilation  and  creation  operators  for  the  slowly  vary¬ 


ing  Stokes  and  anti-Stokes  fields  is  described  by  the  coupled  propagation  equations 


(4.5) 


The  following  analysis  of  the  system  involves  calculation  of  the  expected  values 


of  quantum  field  operators  and  their  combinations.  We  note  that  in  the  Heisenberg 
picture  operators  evolve  and  the  system  is  always  in  its  initial  state,  which  in  our  case 
corresponds  to  no  Stokes  and  anti-Stokes  input  beams  at  the  left  and  right  boundaries. 


4.3  Solution  of  Coupled  Equations 


Due  to  complexity  and  nonlinearity  it  is  not  possible  to  obtain  an  analytic  solution 
for  the  combined  set  of  the  field  propagation  equations  and  Heisenberg-Langevin 


equations.  Nevertheless,  under  the  condition  that  Stokes  and  anti-Stokes  fields  are 
much  weaker  than  the  coupling  and  pump  fields,  and  Stokes  and  anti-Stokes  photon 
densities  are  much  smaller  than  the  atomic  density  N/V  [45],  the  Heisenberg-Langevin 
equations  can  be  linearized  with  Stokes  and  anti-Stokes  aas  fields  as  perturbation 
parameters.  The  linearization  procedure  is  described  in  detail  in  the  Appendix  4.7.3. 

In  order  to  solve  the  set  of  the  linearized  Heisenberg-Langevin  Eq.  (4.35)  and 
coupled  propagation  Eq.  (4.5),  we  first  Fourier  transform  them.  Then,  extracting  the 
solutions  for  and  substituting  them  into  Fourier  transformed  Eq.(4.5), 

we  obtain  the  coupled  equations  for  as(^,a;)  and  —uj)  in  the  form 


da 


(4.6) 
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where  gR{u!),  ras(a;)  are  the  Raman  gain,  EIT  profile  coefficients  respectively,  Ks{ij) 
and  K,as{i^)  are  Stokes  and  anti-Stokes  mode  coupling  coefficients.  In  Eq.  4.6  fcn{z,u) 
are  the  renormalized  Langevin  noise  operators:  fat{z,uj)  =  \/N/c  x  Fcn{z,u)).  The 
sum  is  taken  over  the  relevant  Langevin  noise  operators  {/21, /24, /31, /34}. 

The  general  solution  of  Eq.  (4.6)  can  be  written  as 


where  A/ =  (“/r) 

For  the  following  derivation,  let  us  define  the  coefficients  of  matrix  as 


^1  _  „-ML 

Cl  dJ 


(4.8) 


Due  to  the  linearity  of  Eq.  (4.7),  the  unknown  variables  al{L)  and  aas(O)  of  the 
backward  wave  problem  can  be  written  as  a  linear  combination  of  the  initial  boundary 
values  and  the  noise  terms: 


(4.11) 


The  coefficients  A,B,C  and  D  are  the  functions  of  cj,  whereas  Pai.Qai  ^re  the 
functions  of  uj  and 
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4.4  Characteristics  of  generated  photons 


4.4.1  Stokes  and  Anti-Stokes  photon  generation  rates 

We  evaluate  the  output  generation  rates  of  Stokes  and  anti-Stokes  photons  into  a 
single  transverse  mode,  for  example,  into  a  pair  of  mode-matched  optical  fibers.  The 
generation  rates  at  the  corresponding  boundaries  Zj  are 

(4-12) 

where  subscript  j  denotes  either  Stokes  or  anti-Stokes  photon,  Zs  =  L  and  Zas  =  0. 

Of  importance  are  the  spectral  properties  of  the  generated  photons.  The  power 
spectrum  of  the  output  Stokes  and  anti-Stokes  fields  are  related  to  their  first  order 
coherence  functions  =  (aj(zj^t)aj(zj,t  +  r)}  as 

4-00 

Rj(w)  =  l  j  dTe<^Gf(T)  (4.13) 

— oo 

We  use  the  solutions  for  as(L,  lo)  and  aJs(0,  —lo)  field  operators,  given  by  Eq.  (4.9), 
the  commutation  relations  for  the  input  field  operators  [aj{zj,uj),a}j{zj,  =  L/{2ttc)5{u}+ 
u')  and  Eq.  (4.41).  We  apply  inverse  Fourier  transformation  a{t)  —  J  dwe~'^^a{ijj) 
and  a^(t)  =  f  doje~^‘^^a^—uj)  and  obtain  the  Stokes  and  anti-Stokes  generation  rates 
from  Eq.  (4.12)  in  the  form 


af,ay  Q 

L 


ai,aj  Q 


(4.14) 

(4.15) 


The  integrands  of  Eq.  (4.14)  and  Eq.  (4.15)  are  Stokes  and  anti-Stokes  spectral 
generation  rates  respectively. 

As  seen  from  Eq.  (4.14)  and  Eq.  (4.15),  Stokes  and  anti-Stokes  spectral  generation 
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rates  consist  of  parametric  counts,  characterized  by  the  transfer  functions  C{(jj)  and 
and  of  noise  counts,  that  originate  from  Langevin  noise  fluctuations.  If  the 
pump  is  weak  and  far  detuned,  B{lj)  =  — C(cj),  and  therefore  the  parametric  terms 
in  Stokes  and  anti-Stokes  rates  are  equal.  We  also  note  that  the  contribution  of 
Langevin  noise  fluctuations  to  spectral  generation  rates  can  not  be  neglected  and 
becomes  dominant  at  low  optical  depth  of  the  atomic  sample. 

Using  the  commutator  conservation  relation  for  the  Stokes  held  at  the  right  bound¬ 
ary  z  =  L,  obtained  from  Eq.  (4.9),  the  Stokes  spectral  generation  rate  can  be  ex¬ 
pressed  as 


L 

R,{w)  =  |/1|"  -  1  +  I  (4.16) 

ai,aj  Q 

When  the  pump  is  far  detuned  from  the  atomic  transition  and  Stokes  photon  losses 
are  small,  the  contribution  of  Langevin  noise  fluctuation  in  Eq.  (4.16)  is  negligible. 
Thus,  the  Stokes  spectral  generation  rate  can  be  written  as 

(4.17) 

Eq.  (4.17)  can  be  interpreted  in  terms  of  the  quantum  theory  of  linear  amplifica¬ 
tion  [47]:  in  the  absence  of  the  Stokes  input  beam  and  losses  for  the  Stokes  photon,  the 
Stokes  generation  rate  is  just  the  additive  noise  caused  by  the  amplification  process, 
that  is  characterized  by  the  gain  coefficient  A{uj). 

4.4.2  Two-photon  intensity  correlation  function  and  bipho¬ 
ton  wave-function 

In  order  to  address  another  important  issue  -  the  time  correlation  properties  of  the 
generated  photons,  we  calculate  the  Glauber  two-photon  correlation  function  of  time 
delay  r  between  Stokes  and  anti-Stokes  photons: 


<^i-L(T)  =  ^  +  T)«as(0,  t  +  T)as{L,  t)) 


(4.18) 
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Using  the  Stokes  and  anti-Stokes  operators  at  the  boundaries  derived  earlier 
[Eq.  (4.9)]  and  the  commutator  relations  for  the  input  fields,  we  obtain  the  inten¬ 
sity  correlation  function  as 


G 


(2) 

s—as 


("t) 


X  -wi)aL(0-  -W2)aas(0, 0)3)05(1, 0)4)) 


(4.19) 


The  intensity  correlation  function  contains  the  fourth-order  Langevin  noise  cor¬ 
relations.  According  to  Gaussian  moment  theorem  [48,  49]  they  can  be  decomposed 
preserving  the  order  into  the  sum  of  the  products  of  second  order  Langevin  noise 
correlations.  As  a  result  the  intensity  correlation  function  can  be  simplified  to 


=  (4s(0>  ^  +  'r)aa5(0,  t  -b  r))  {al{L,  t)as{L,  t))  -b  |  {aas{0,  t  +  T)as{LH))  f 
=  G(i)(0)  X  GW(0)  +  \^s-as{r)f  (4.20) 

The  first  term  in  Eq.  (4.20)  represents  flat  uncorrelated  background,  the  second 
term,  expressed  through  ^s-as{^)  function,  describes  the  correlation  part.  ^s-as{^) 
is  equal  to 


^s-as{r)  =  {aas{0,t  +  r)as{L,t)) 

L 


L 

2ttc 


f  dwe^^(BD*  +Y^  f 

ai,aj  Q 


(4.21) 


We  note  that  e  represents  a  two  photon  wavefunction  or  a 

biphoton  wavefunction  on  the  condition  that  peak  value  of  normalized  g^^las  >>  1j 
where  gflasiT)  =  Gflas{r) / {G^^\o)  x  ^1^(0)). 

Similarly,  we  can  define  4>as-s(T)  =  {as{L,t)aas{0,t  +  r)).  This  function  can  be 
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obtained  in  the  form 


$as-s(T)  = 


(4.22) 


L 


L 

27rc 


For  a  wide  range  of  input  parameters  we  numerically  verify  that  4>a5_s(T)  = 


For  the  case  where  Qp/  AuJu  «  I  and  AuJuhu  »  1,  we  numerically  verify  that 
the  contribution  of  Langevin  noise  fluctuations  to  4>a5_5(T)  is  negligible,  therefore 
Eq.  (4.22)  can  be  simplified  to 


(4.23) 


4.5  Ideal  Spontaneous  Parametric  Down  Converter 


Before  we  proceed  to  the  discussion  of  the  interesting  cases  of  the  atomic  correlation 
functions  and  photon  spectral  densities,  we  want  to  make  an  analogy  to  the  well- 


known  parametric  down  converter  in  crystals  [40].  We  consider  the  ideal  model  - 


non-degenerate  parametric  down  converter  in  which  a  generated  signal  photon  has  a 
very  slow  group  velocity  Vg  as  compared  to  an  idler  photon.  We  assume  that  both  idler 
and  signal  photon  escape  SPDC  without  losses,  therefore  the  Langevin  noise  terms  in 
Eq.  (4.6)  can  be  neglected.  In  crystals  the  coupling  coefficient  can  be  approximated  as 
a  constant  over  the  broad  spectral  range  K5(a;)  =  ^^5(0;)  =  k.  We  also  neglect  Raman 
gain  and  approximate  EIT  profile  as  r(a;)  =  —iu)jVg.  Under  these  assumptions 

the  signal  and  idler  photons  have  identical  spectral  characteristics  and  rates.  Using 
the  Eq.  (4.14),  we  obtain  the  photon  spectral  density  R{u))  and  spectrally  integrated 
generation  rate  R  =  l/(27r)  f  dujR(uj)  in  the  form: 
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R(u;)  =  iKpL^sinc^ 

R  =  Vg\K\‘^L. 


(4.24) 


As  seen  from  Eq.  (4.24),  the  spectral  bandwidth  of  the  SPDC  in  crystals  is  limited 
by  Au;  ~  27rV^/L  due  to  the  phase  mismatch  of  the  off-centered  counts  resulting  from 
the  group  delay  in  the  signal  channel.  In  principle,  the  bandwidth  can  be  made 
very  small  by  making  the  group  delay  Tg  =  LjVg  large.  For  the  atomic  system  k  is 
proportional  to  the  atomic  density  Af  =  NjV ,  therefore  the  spectral  density  of  the 
generated  photons  scales  as  {AfVf.  Taking  into  account  that  the  EIT  induced  group 
velocity  of  a  wave  is  Vg  =  n^/(27i3A/’cri3),  where  =  Pi3C^i3/(ceo^7i3)  is  the  atomic 
cross  section  of  the  |1)  — >  |3)  transition,  the  total  count  rate  scales  linearly  with  AfL. 


Figure  4.2:  Normalized  signal-idler  intensity  correlation  function  for  the  ideal  SPDC 
with  the  condition  of  a  large  group  delay  L/Vg  for  a  signal  photon.  R  is  the  paired 
photon  generation  rate. 

Fig.  4.2  shows  the  normalized  signal  idler  intensity  correlation  function.  Since 
the  probability  of  emitting  a  photon  pair  is  uniformly  distributed  along  the  crystal 
of  length  L  and  the  signal  photon  has  a  group  delay  relative  to  the  idler,  the  5^i\(r) 
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is  expected  to  be  the  shifted  rectangle  with  the  width  equal  to  the  group  delay  Tg  = 
L/Vg.  By  varying  the  group  delay  we  can  control  the  width  of  the  waveform.  The 
peak  height  of  is  equal  to  1/(kL)^  and  can  be  interpreted  as  a  duty  ratio:  = 

l/{RTg).  In  the  regime  of  very  small  parametric  gain  kL  «  1,  which  can  be  achieved, 
for  example,  by  lowering  the  intensity  of  the  pump,  the  atomic  system  can  produce 
highly  correlated  photon  pairs  >>  1. 

4.6  EIT  based  Paired  Photon  Generator 

Now  we  turn  from  the  discussion  of  the  ideal  SPDC  to  the  discussion  of  the  EIT 
based  paired  photon  generator.  In  the  EIT  based  generator,  coupling  between  Stokes 
and  anti-Stokes  modes  is  bandwidth  limited.  Moreover,  the  generated  Stokes  photon 
undergoes  Raman  gain  whereas  a  paired  anti-Stokes  photon  propagates  slowly 

and  undergoes  absorbtion  at  the  poles  of  EIT  profile  a;  =  ±£1^2. 

We  first  obtain  the  coefficients  of  Eq.  (4.6)  with  the  assumption  that  the  pump  is 
weak  and  far  detuned  from  the  |1)  — |4)  transition  and  Afc  =  {kp+kc  —  ks  —  kas)'Z  =  0: 


CHAPTER  4.  THEORY  OF  EIT  BASED  PAIRED  PHOTON  GENERATION  42 


r  = 

A  aA  — 


2iAfa^n{u)  +  ^712) 


9r 


G{u) 

^  \  +  i7i3) 


\2NojIJ 


(  % 

Kg  =  Kas  =  ~  ^  ^  ■ 

\  2NiJ 
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G{u) 
iAfa^n 
G{u)) 


(4.25a) 

(4.25b) 

(4.25c) 
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Dc\/ Afa^iz 

2v^Aa;i4^  ’ 


(4.25d) 

(4.25e) 

(4.25f) 

(4.25g) 


^21  = 


\/2J7c\/A/"(77i3 


tas  _  _ 

S24  “ 

Ql  = 

^34  = 


G(a;)  ’ 

Op  \  Oc\/ A/'(T713 


AtJi4  y  \/2G'(tj) 

2\/2(a;  +  i7i2)\/A/o^ 

G(a;) 

^  \/2(a;  +  i7i2)\/A/'tT7i3 
Aa;i4  y  G(a;) 


(4.25h) 


(4.25i) 


(4.25j) 

(4.25k) 


where  A/"  is  an  atom  density,  G{u})  =  If^d^  —  4(a;  +  z7i2)(a;  +  i7i3),  a  is  the  absorption 
cross  section  for  all  allowed  transitions:  a  —  uu  —  ^24  =  <723  =  o'i3- 
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4.6.1  Ground  state  approximation 

The  approximation  that  the  atomic  population  remains  in  the  ground  state  an  =  1 
gives  two  significant  diffusion  coefficients  Pi2,2i  =  2712  and  ^13,31  =  r3  correspond¬ 
ing  to  ^21,^31  Langevin  noise  operators.  Due  to  small  atomic  population  in  the 
excited  states,  the  rest  of  the  diffusion  coefficients  are  approximated  as  zeros  and  the 
corresponding  Langevin  noise  operators  are  neglected  [36]. 

We  numerically  examine  the  emission  rates  and  intensity  correlation  function  for 
the  EIT  based  paired  photon  emitter.  We  take  other  parameters  similar  to  those 
of  a  Rb  MOT:  atom  density  J\f  =  10^^  atoms  per  cm^,  atomic  cross  sections  a  = 
cri3  =  ai4  =  a24  =  10^^  cm^  and  dephasing  rates  equal  to  one  half  of  the  Einstein  A 
coefficient,  i.e.,  713  =  714  =  724  =  723  =  L79  x  10^  radians.  We  choose  the  strength 
and  the  detuning  of  the  pump  laser  from  the  |1)  ^  |4)  transition  as  AuJu  =  24713 
and  Qp/Auju  =  0.1. 


Figure  4.3: 


Transmission  and  phase  mismatch  as  functions  of  detuning  uj 


Fig.  4.3  shows  the  EIT  transmission  profile  and  phase  mismatch  as  a  function  of 
the  detuning  of  the  anti-Stokes  frequency  u.  We  take  Qc  =  6713,  AfaL  =  11  and 
7i2  =  0.  Fig.  4.4  shows  the  profiles  of  the  coupling  coefficient  |/^(ct;)|  and  the  Raman 
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Figure  4.4:  Coupling  constant  and  Raman  gain  as  functions  of  detuning  uj. 

gain  coefficient  Re{gfi{uj))  for  the  same  parameters  as  in  Fig.  4.3. 

EIT  based  spontaneous  emitters  allow  a  variable  emission  bandwidth.  The  band¬ 
width  and  profile  of  the  spontaneous  emission  rates  are  controlled  by  the  strength  of 
the  coupling  laser  Qc  ^nd  the  optical  depth  of  the  atomic  sample  AfaL.,By  reducing 
the  strength  of  the  coupling  laser,  the  emission  linewidth  can  be  made  much  smaller 
than  the  natural  linewidth,  with  the  minimum  width  ultimately  limited  by  the  de¬ 
phasing  rate  712  of  the  |1)  — >  |2)  transition.  At  small  optical  depth  AfaL  <  1,  the 
Raman  gain  coefficient  and  the  coupling  coefficients  Ks{uj)^  ^as(^)  determine 

the  emission  spectrum.  At  high  optical  depth  AfaL  >  1,  the  EIT  transmission  win¬ 
dow  and  the  phase  mismatch,  introduced  by  a  large  group  delay  in  the  anti-Stokes 
channel,  affect  the  spontaneous  emission  spectrum. 

In  Fig.  4.5,  Fig.  4.6  and  Fig.  4.7  we  show  the  variations  of  the  coincidence  count 
rate  Rc{t)  in  a  1  ns  bin  and  the  corresponding  Stokes  power  spectral  density  depend¬ 
ing  on  the  optical  depth.  With  a  bin  size  AT  =  1  ns  much  smaller  than  the  correlation 
time,  the  coincidence  count  rate  is  obtained  from  the  intensity  correlation  function 
as  Rc{t)  =  AT{c/LyG^as-s{'^)'  Compared  to  the  ideal  case,  described  earlier,  the 
intensity  correlation  function  and  the  emission  spectrum  for  the  EIT  based  atomic 
SPDC  show  some  interesting  features.  The  shape  of  the  intensity  correlation  function 


CHAPTER  4.  THEORY  OF  EIT  BASED  PAIRED  PHOTON  GENERATION  45 


Figure  4.5:  The  oscillatory  regime:  (a)  Coincidence  count  rate  in  a  1  ns  bin  and  (b) 
Stokes  spectral  generation  rate.  AfoL  =  0.3,  Dc  =  6713  and  712  =  0. 


Figure  4.6:  The  group  delay  regime:  (a)  Coincidence  count  rate  in  a  1  ns  bin  and  (b) 
Stokes  spectral  generation  rate.  MaL  =  20,  Qc  =  6713  and  712  =  0. 
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Figure  4.7:  The  group  delay  regime:  (a)  Coincidence  count  rate  in  a  1  ns  bin  and  (b) 
Stokes  spectral  generation  rate.  AfcrL  =  200,  Oc  —  6713,  Dp  =  0.8713  and  712  =  0. 
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Figure  4.8:  Intensity  correlation  function  with  the  correlation  time  smaller  than  the 
spontaneous  decay  time,  Tr  <  Tg  <  1/(2713).  J\faL  =  200,  flc  =  66713,  712  =  0  and 
1/(2713)  =  25  ns. 


CHAPTER  4.  THEORY  OF  EIT  BASED  PAIRED  PHOTON  GENERATION  47 


and  the  emission  profile  depend  on  the  relation  of  three  characteristic  times.  The  first 
is  the  inverse  Rabi  frequency  of  the  coupling  laser  =  2;r /  —  713 ,  the  second  is 

the  group  delay  between  Stokes  and  anti-Stokes  photons  Tg  =  L/Vg  =  2'yi3NfcrL/Q.‘^ 
and  the  third  is  a  minimal  pulse  length  required  to  pass  through  EIT  medium  [50] 

rp(min)  =  81n(2)7i3\/A/'(TL/n2. 

When  the  EIT  effect  is  small,  which  occurs,  for  example,  at  low  optical  depth,  the 
atomic  system  behaves  like  a  single  atom  [51].  In  such  a  regime  the  intensity  corre¬ 
lation  function  reveals  the  damped  Rabi  oscillations  (Fig.  4.5).  The  oscillations  ob¬ 
served  in  the  intensity  correlation  function  have  the  time  period  of  =  27r/  ~  7i3 

and  occur  on  the  condition  that  >  Tg,  Tp^rnin)  and  the  coupling  laser  is  strong  enough 
to  force  the  oscillations  to  overcome  damping  ilc  >  7i3-  Once  the  metastable  state  |2) 
is  excited  by  the  Raman  process  |1)  |4)  |2),  the  probability  amplitude  between 

1 2)  and  |3)  oscillates  due  to  the  strong  interaction  of  the  atoms  with  the  resonant 
coupling  beam. 

In  Fig.  4.6  we  show  the  intensity  correlation  function  and  Stokes  emission  spectrum 
in  the  group  delay  regime,  where  Tg  >  Tr  and  Tg  >  Tp^min)-  In  this  regime  the  width 
of  the  intensity  correlation  function  is  approximated  by  Tg.  Moreover,  the  frequency 
range  over  which  the  spontaneous  generation  occurs  is  filtered  by  the  EIT  window 
and  mostly  controlled  by  the  phase-matching  in  the  presence  of  large  group  delay 
in  the  anti-Stokes  channel.  A  sufficiently  wide  EIT  window  Tg  >  Tp^min)  requires 
high  optical  depth  AfaL  >  10.  In  the  presence  of  the  non- zero  dephasing  712,  the 
maximum  group  delay  and  therefore  the  maximum  width  of  the  correlation  function 
is  limited  by  1/712- 

By  increasing  the  optical  depth  of  the  atomic  sample,  the  EIT  window  can  be  made 
substantially  larger  than  the  emission  bandwidth.  A  large  EIT  window  might  be  very 
useful  for  such  applications  as  bi-photon  waveform  control  and  shaping.  Fig.  4.7 
shows  the  intensity  correlation  function  and  the  Stokes  emission  spectrum  under  such 
a  condition,  where  AfcrL  =  200.  The  tail  of  the  correlation  function  decays  on  the 
time  scale  of  Tp(mm)- 

We  also  note  that  at  high  optical  depth  it  becomes  possible  to  achieve  the  Stokes- 
anti-Stokes  correlation  function  with  the  width  shorter  than  the  spontaneous  decay 
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time  1/(2713).  Fig.  4.8  shows  such  a  correlation  function  that  is  obtained  for  NaL  = 
200  and  Qc  =  66713,  where  Tr  <  Tg  <  1/(2713). 

We  turn  next  to  the  sharp  peak  at  the  leading  edge  of  the  correlation  function  of 
Fig.  4.7.  This  sharp  peak  is  Sommerfeld-Brillouin  precursor  [52],  It  has  been  observed 
in  our  theoretical  plots  for  a  long  time.  But  only  recently,  when  it  was  experimentally 
measured,  we  got  interested  in  its  nature.  We  thank  Daniel  J.  Gauthier  who  suggested 
that  this  peak  is  the  result  of  simultaneously  generated  Stokes  and  anti-Stokes  photons 
that  travel  at  nearly  the  speed  of  light  in  vacuum  and  arrive  near-simultaneously  at 
the  photodetectors.  He  also  pointed  out  to  the  similarity  of  this  peak  to  precursors 
which  have  been  extensively  studied  [52,  53,  54].  A  precursor  has  an  approximate 
width  that  is  equal  to  the  opacity  width  of  the  atomic  transition  in  the  optically  thick 
medium.  Similar  to  precursors,  the  sharp  leading  edge  peak  has  a  opacity  width  of 
the  EIT  profile  [55]. 

4.6.2  No  ground  state  approximation 

If  the  pump  is  weak  and  far  detuned  and  therefore  most  of  the  atomic  population  is 
in  the  ground  state,  we  verify  that  the  ground  state  approximation  gives  a  correct 
prediction  for  the  Bi-photon  function  and  the  Stokes  generation  rate.  Nevertheless 
it  does  not  properly  account  for  an  atom  return  to  the  ground  state  |1).  Ideally  one 
would  expect  the  Stokes  and  anti-Stokes  rates  to  be  equal,  since  an  atom,  making  a 
complete  cycle  on  the  energy  level  diagram  (Fig.  5.7),  returns  to  the  ground  state. 
Even  at  zero  dephasiiig  rate  712  =  0  of  the  |1)  ^  |2)  transition,  the  ground  state 
approximation  predicts  the  anti-Stokes  generation  rate  to  be  smaller  than  the  Stokes 
generation  rate.  For  example,  Ras/Rs  =  0.65  for  AfaL  =  10,  Qp/Auji^  =  0.1  and 
Oc  =  6713. 

In  order  to  treat  properly  an  atom’s  return  to  the  ground  state  we  will  retain  in 
Eq.  (4.6)  all  four  Langevin  noise  operators  {/21, /24? /31? /34}  and  take  into  account 
small  incoherent  population  in  excited  states,  resulting  from  the  steady  state  solutions 
of  Eqs.  (4.27a)-(4.27f).  With  these  inclusions  the  solution  of  Eq.  (4.6)  predicts  the 
Stokes  and  anti-Stokes  spectral  generation  rates  to  be  equal  at  a  zero  dephasing  rate 
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7i2  =  0  [Fig.  4.9(a)].  A  non-zero  dephasing  rate  712  ^  0  reduces  EIT  and  therefore 
introduces  additional  lossas  for  anti-Stokes  photons.  As  a  result  the  output  Stokes 
rate  exceeds  the  anti-Stokes  rate.  For  example,  for  AffjL  =  10,  Qp/Auu  =  0.1, 
=  6713  and  7i2  =  O.6713  we  obtain  Ras/Rs  ~  0.8.  The  corresponding  Stokes  and 
anti-Stokes  spectral  generation  rates  are  shown  in  Fig.  4.9(b).  We  believe  that  the 
“missing”  anti-Stokes  photons  are  absorbed  and  are  then  reemitted  in  a  solid  angle 
of  47r.  The  atomic  sample  is  assumed  to  be  optically  thin  in  the  radial  direction. 


Figure  4.9:  Stokes  (solid  curve)  and  anti-Stokes  (dash  curve)  spectral  generation  rates 
at  (a)  zero  dephasing  712  =  0  and  (b)  non  zero  dephasing  712  =  O.6713.  AfaL  =  10, 
=  6713  Dp  =  2.4713,  Au;i4  =  24713 


Compared  to  the  ideal  SPDC,  where  each  generated  signal  photon  has  its  paired 
idler  photon,  the  real  atomic  system  has  uncorrelated  noise  counts  in  both  Stokes 
and  anti-Stokes  channels  that  result  from  Langevin  noise  fluctuations.  In  Fig.  4.10 
we  examine  the  dependence  of  the  Stokes  (anti-Stokes)  and  paired  count  rates  on 
the  optical  depth  AfaL.  The  parameters  for  the  curves  are  712  =  0,  Qc  =  5713, 
Au;i4  =  24713  and  Dp/AuJu  =  0.1.  The  paired  count  rate  {Rp)  is  defined  as  the  area 
under  the  Stokes-anti-Stokes  coincidence  count  rate  function  minus  the  area  under 
the  uncorrelated  background.  One  may  show  that  Rp  ^  V(27r)  /  duJ  lAC*]^.  At  small 
optical  depth  the  paired  rate  scales  quadratically  with  the  optical  depth  and  is  much 
smaller  than  the  Stokes  rate.  At  high  optical  depth  the  paired  rate  varies  linearly 
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Optical  Depth 


Figure  4.10:  Stokes,  anti-Stokes  (solid  curve)  and  Paired  (dashed  curve)  photon  gen¬ 
eration  rates  as  a  function  of  the  optical  depth.  Qc  =  ^713,  712  =  0,  Acc;i4  =  24713  and 
Dp/Auji4  =  0.1.  At  optical  depth  of  100,  the  paired  rate  reaches  90%  of  the  Stokes 
rate. 

with  J\fcrL  and  converges  logarithmically  to  the  Stokes  emission  rate. 

In  this  chapter  we  describe  the  theory  of  paired  photon  generation  in  double- 
A  atomic  system.  With  low  parametric  gain  and  high  optical  depth  we  show  that 
the  system  can  produce  highly  correlated  photon  pairs.  The  shape  of  the  intensity 
correlation  function  and  the  emission  bandwidth  depend  on  the  coupling  laser  Rabi 
frequency  and  the  optical  depth  of  the  atomic  sample.  Compared  to  the  ideal  SPDC, 
paired  photon  generation  in  the  double-A  atomic  system  is  affected  by  Raman  gain  in 
the  Stokes  channel  and  EIT  in  the  anti-Stokes  channel.  EIT,  through  the  absorption 
at  the  poles,  cuts  the  emission  bandwidth.  In  order  to  enter  a  regime  where  the  EIT 
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window  is  sufficiently  large  and  therefore  the  emission  bandwidth  is  controlled  to  a 
large  extent  by  the  phase-matching  process  in  the  presence  of  large  group  delay,  the 
optical  depth  of  the  atomic  sample  has  to  be  large  AfcrL  >  10.  High  optical  depth 
substantially  reduces  the  influence  of  Langevin  noise  fluctuations  on  paired  photon 
generation  so  that  the  Stokes  and  anti-Stokes  photons  are  generated  mostly  in  pairs. 
We  therefore  suggest  the  use  of  a  cigar  shaped  atomic  cloud  with  high  optical  depth 
in  the  longitudinal  direction. 

4.7  Mathematical  Apparatus 

4.7.1  Collective  slowly  varying  atomic  operators 

To  describe  the  quantum  properties  of  the  atomic  system  we  use  the  collective  slowly 
varying  atomic  operators  [56,  57,  58,  46,  59]  d'jk(z,t),  deflned  as 

i)  =  >  (4-26) 

^  i&Nz 

where  the  averaging  is  done  over  each  atom  z  in  a  small  interval  Az  that  contains 
large  number  of  atoms  Nz  »  1.  The  slowly  varying  variables  are  assumed  to  stay 
unchanged  over  Az.  —  coi  +  AuJu,  u^2  =  W4  —  ozs  +  AuJu,  Uu  =  —  i^4i, 

1^24  =  — i^42,  the  rest  of  Ujk  —  ujj  —  oJk-  hi  —  hs  ■  z,  ^42  —  kg-z,  kn  =  kp-z,  k^^  =  kc-z 
are  the  projections  of  the  anti-Stokes,  Stokes,  pump  and  coupling  /c-vectors  on  the  2 
axis,  k^3  =  A;4i  —  k^i,  k2i  =  —k2i  +  k4^i,  the  rest  of  kjk  =  —k^j. 

4.7.2  Heisenberg-Langevin  equations 

The  full  set  of  Heisenberg-Langevin  equations  for  the  four  state  system  consists  of  16 
equations.  Here,  we  show  explicitly  10  of  them,  the  other  6  for  the  adjoint  off-diagonal 
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atomic  operators  are  not  shown. 


^5-44  =  F44  —  r4a44 
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(4.27a) 


(4.27b) 


(4.27c) 


(4.27d) 


(4.27e) 


(4.27f) 
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d 

— a'34  =  F34  —  (734  —  iNujii)  a34 
at 

Al  (  ”^<^31  +  9s<1s(^32  ~  dasO'as^U - ^*^24 

a,  ~ 
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—d'2A  =  F^i  —  (724  ~  iAa;i4)  a24 
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(4.28a) 


(4.28b) 


(4.28c) 


(4.28d) 


Here,  for  simplicity,  we  assume  that  Ak  =  {kp  +  kc  —  ks~  kas)  ■  z  =  0.  In  Eq.  (4.27a)- 
Eq.  (4.28<i)  Ej  is  the  total  decay  rate  from  state  |i),  Ej-j  is  the  decay  rate  from  state 
|i)  to  state  |j)  and  7,j  is  the  dephasing  rate  between  state  |i)  and  state  |j).  The 
dephasing  rates  for  the  double-A  system  in  the  absence  of  the  collisional  dephasing 
can  be  obtained  from  total  decay  rates  E3  and  E4  from  state  |3)  and  |4)  to  two  ground 
states  |1)  and  |2)  as 


731 


732  =  Y 


E4 

741  =  742  =  Y 

E3  +  E4 


743  = 


2 


(4.29) 

(4.30) 

(4.31) 
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4.7.3  Linearization  procedure 

In  zeroth  order  perturbation  expansion,  in  which  Og  and  has  go  to  zero,  the  Heisenberg- 
Langevin  equations  for  hii,ai4,a22,h23,d’32,h33,a4i,a44  atomic  operators  are  decou¬ 
pled.  Under  the  assumption  that  pump  and  coupling  beams  propagate  without  de¬ 
pletion,  we  obtain  the  steady  state  solution  for  Eqs.  (4.27a)-(4.27f)  in  form 

=  {^%)  +  ^  (mnFmn  (4.32) 

With  the  definition  of  the  denominator  as 

T  =  Tsi  (r^  -p  +  2  |f7p|") 

+  r42(rl  +  2|f7/)  (4.33) 

the  steady  state  expectation  values  for  the  zeroth  order  atomic  operators  are  equal 
to 


/.o\  r3i|f7,|"(r2  +  4Aa;?4  +  |f7pp) 

V  11/  2^ 

(4.34a) 

r„(ri  +  |n/)|!y" 

W22/  -  j. 

(4.34b) 

/-O  \  _  ^42 

\^33/  rp 

(4.34c) 

/-O  \  _  ^31 
\^44/  p 

(4.34d) 

/-o  \  _  Tsi  (2Aa;i4  —  fr4)  iricl^  ^^p 

Wl4/  p 

(4.34e) 

/-o  \ _ ^r3r42ric  |f^p| 

\  ^23  /  p  ’ 

(4.34f) 

In  the  first  order  expansion,  we  substitute  the  zeroth  order  solution  for  the  atomic 
operators  Eq.  (4.32)  into  the  remaining  Heisenberg- Lange vin  equations  for  0^21,^24, 
^31  ^  ^34  and  their  adjoint.  Neglecting  higher  order  terms  like  emnFmn^s  or  emn^mnSL 
we  obtain  the  linearized  equations.  We  note  that  the  linearized  Eqs.(4.28a)-(4.28d) 
for  <521,^24,5-31,(734  represent  an  independent  set  of  equations  and  can  be  decoupled. 
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For  clarity  we  write  them  in  vector  form 

— =  yAo'i  +  +  Fi  (4.35) 

where  ai  =  {^21,^24,^31,^34},  Fi  =  {F21, ^4,  Ai,  A4},  a  =  {as,aL})  matrix  A 
depends  on  the  dephasing  rates  7^^  and  pump  and  coupling  laser  rabi  frequencies 
fip,  Clci  matrix  A4  depends  on  the  zeroth  order  solution  for  the  atomic  operators 

{^%)- 

4.7.4  Langevin  noise  operators  and  their  diffusion  coefficients 

By  analogy  with  the  collective  slowly  varying  atomic  operators  ajk{z^  t)^  the  collective 
Langevin  noise  operators  are  defined  as 

^  Zi^Nz 

We  assume  that  a  Langevin  noise  operator  for  a  single  atom  is  5-correlated  so  that 

{FSit)F*!lit'))  =  (4-37) 

where  (...)  denotes  the  average  over  the  reservoir,  is  the  atomic  diffusion 

coefficient  for  an  ith  atom. 

Now  we  consider  the  second  order  correlations  for  the  collective  Langevin  noise 
operators 

{Ft,{t.z)FMt'.^'))  =  4  E  (4.38) 

^  Zi&N, 

Introducing  the  average  atomic  diffusion  coefficient 

Djk,i'k'{l,z)  =  — 

^  Zi^N, 


(4.39) 
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the  noise  correlations  can  be  expressed  as 

z)Fj>k’{t',  z'))  =  z)5{t  -  t')5{z  -  z')  (4.40) 

In  case  Djkj'k'{Iy  is  independent  of  t,  in  frequency  domain  the  noise  correlations 
are 

{Fjk{uj,  z)Fj'k'{uj',  z'))  =  i^Rjkj'k'^iu)  +  uj')5{z  -  z')  (4.41) 

The  diffusion  coefficients  Djkj'k'  can  be  obtained  from  the  Heisenberg-Langevin 
equations  (4.27a)-(4.28d)  using  the  generalized  fluctuation-dissipation  theorem  [60, 
61].  Here  we  show  the  diffusion  coefficients  for  the  Langevin  noise  operators  of  interest 
T21,T24,T3i,F34  and  their  adjoint  ^12,^42,^13,^43 
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(4.43) 

where  denotes  {21,24,31,34}  subspace  for  the  atomic  operators,  a]  denotes 


(12,42, 13,43}  subspace  for  the  adjoint  atomic  operators. 


